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Exact spectrum of the Laplacian on a 
domain in the Sierpinski gasket 

HUA QIU 

Abstract. For a certain domain Q in the Sierpinski gasket SQ whose boundary is a hne 
segment, a complete description of the eigenvalues of the Laplacian under the Dirichlet 
pH , and Neumann boundary conditions is presented. The method developed in this paper is a 

^ ■ weak version of the spectral decimation method due to Fukushima and Shima, since for a 

■ lot of "bad" eigenvalues the spectral decimation method can not be used directly. We also 

prove an analogue of Weyl's classical result on the eigenvalue asymptotics of the eigenvalue 
< : counting function The ratio p^(x)/xi°s3/iog5 bounded but non- convergent as 

X — oo. Moreover, we explain that the asymptotic expansion of admits a second 

term of the order log 2/ log 5, that becomes apparent from the experimental data. This is 
very analogous to the conjectures of Weyl and Berry. 

Keywords. Sierpinski gasket, Laplacian, eigenvalues, spectral decimation, analysis on 
fractals. 
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1 Introduction 

^ , The study of the Laplacian on fractals was originated by S. Kusuoka [21] and S. Goldstein 
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[TT] . They independently constructed the Laplacian as the generator of a diffusion process 
on the Sierpinski gasket SQ. Later an analytic approach was developed by J. Kigami, 
who constructed the Laplacian both as a renormalized limit of difference operators and a 
weak formulation using the theory of Dirichlet forms [15] . 

We are particularly interested in the eigenvalues of this Laplacian. In the case of the 
Sierpinski gasket. Physicists R. Rammal and G. Toulouse [27] found that an appropriate 
choice of a series of eigenvalues of successive difference operators produces an orbit of 
the dynamical system related to a quadratic polynomial. This is the phenomenon which 
M. Fukushima and T. Shima [lOl |30] described from the mathematical point of view, 
by saying that SQ admits spectral decimation with respect to a quadratic polynomial. 
Furthermore, they found all the eigenvalues of the Laplacian on the Sierpinski gasket by 
tracking back the orbits. Later the theory of the Laplacian was developed for nested 



This research was supported by the National Science Foundation of China, Grant 10901081. 



1 



fractals and p.c.f. self-similar sets by T. Lindstr0m [23] and J. Kigami [16] by introducing 
the notion of harmonic structure. Every p.c.f. self-similar set is approximated by an 
increasing sequence of finite graphs and the harmonic structure determines a sequence of 
difference operators on the successive graphs, which converges to the Laplacian. Then 
some generalizations of the spectral decimation to a class of p.c.f. self-similar sets were 
developed by T. Shima [31], L. Malozemov and A. Teplyaev [2l], in which some strong 
symmetry conditions are supposed to be satisfied to ensure the spectral decimation applies 
to the corresponding graph sequences. Under such strong symmetry conditions, the spec- 
trum of the Laplacian can also be described in terms of the iteration of a rational function. 
Recently, the spectra of the Laplacian operators on some other fractals have been ana- 
lyzed either numerically [1] or using the spectral decimation method [H [H [371 EH] by R. 
S. Strichartz (with co-authors), D. Zhou and A. Teplyaev. In all the references mentioned 
above, spectral decimation plays a key role in the theoretical study of the spectra of the 
Laplacian operators. 

In this paper, we are mainly concerned with eigenvalue problems for a domain in the 
Sierpinski gasket. Since analysis on fractals has been made possible by the definition 
of Laplacian, it is natural to explore the properties of these fractal Laplacians that are 
natural analogs of results that are known for the usual Laplacian. However, since not 
much is known about the fractal Laplacian, we can only scratch the surface in attempting 
the generalization to fractal Laplacian. 

For simplicity, here we specifically focus on the Sierpinski gasket SQ. Recall that SQ 
is the attractor of the iterated function system {Fq, Fi, F2} where FiX = |(x + qi) where 
90) 9i) Q2 are the vertices of an equilateral triangle in the plane, 

2 

sg = [jFi{sg). 

i=0 

Let A denote the Laplacian on SQ defined by Kigami. In his theory the boundary of 
consists of the three points qo,qi, q2 and the space of harmonic functions (solutions of 
An = 0) is three dimensional, with u determined explicitly by its boundary values u{qi). 
(Note that this boundary is not a topological boundary.) Thus this theory is more closely 
related to the theory of linear functions on the unit interval than to harmonic functions 
on the disk. To get a richer theory we should take an open set in iS^ and restrict 
the Laplacian on SQ to functions defined on Q. Thus harmonic functions on Q are the 
solutions of An = 0. 

Here we particularly focus on the certain domain fi^. which is a triangle obtained by 
cutting SQ with a horizontal line at any vertical height x (0 < x < 1 if we suppose that 
the height of SQ is equal to 1.) below the top vertex go- See Fig. 1.1. An important 
motivation for studying this kind of domains is that they are the simplest examples which 
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Fig. 1.1. Qx and fli. 



could serve as a testing ground for questions and conjectures on analysis on more general 
fractal domains with fractal boundaries. These domains were first introduced by R. S. 
Strichartz in [32] and later studied by J. Owen and R. S. Strichartz in [25], where they 
gave an explicit analog of the Poisson integral formula to recover a harmonic function u 
on Vtx from its boundary values. It is also natural to calculate an explicit Green's function 
for the Laplacian on VL^. This was studied by Z. Guo, R. Kogan and R. S. Strichartz in 
[12] which is completely similar to the construction of the Green's function on SQ given 
by Kigami in [T5l 116] . For some other analytic topics related to this kind of domains, see 

[ni [HI na [20]. 

In the present paper, we study the spectral properties of the Laplacian on Vt^ which is 
an open problem posed in [25]. For the simplicity of description, we mainly concentrate 
our attention to a particular domain Vti (We drop the subscript 1 on in all that follows 
without causing any confusion.) which is the complement of go U where L is the line 
segment joining qi and q2 (in this case dVt = qqU L). We give a complete description of 
the Dirichlet and Neumann spectra of the Laplaician on Q. 

Unfortunately in our context, for a number of "bad" eigenvalues (whose associated 
eigenfunctions have supports touching the bottom boundary line L) spectral decimation 
can not be used directly, which makes things more complicated. However, by choosing 
a sequence of appropriate graph approximations, we describe a phenomenon on those 
eigenvalues called weak spectral decimation which approximates to spectral decimation 
when the levels of the successive graphs go to infinity. And we use this weak spectral dec- 
imation to replace the role of spectral decimation in the original Fukushima and Shima's 
work [To]. Actually, similar to the standard case, weak spectral decimation can also pro- 
duce a "weak" orbit related to the same quadratic polynomial by an appropriate series 
of eigenvalues of successive difference operators on graph approximations. We can then 
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trace back those "weak" orbits to capture all the "bad" eigenvalues. More precisely, we 
classify the eigenvalues of A on into three types, which we call the localized eigenvalues, 
primitive eigenvalues and miniaturized eigenvalues. The localized eigenf unctions associ- 
ated to localized eigenvalues on Q are just a subspace of the localized eigenfunctions on 
SQ whose supports are disjoint from L. This kind of eigenvalues can be dealt with in a 
completely similar way to the SQ case, for which the spectral decimation can apply. The 
primitive and miniaturized eigenvalues are the so-called "bad" eigenvalues. They are the 
eigenvalues need to be paid particular attention to. 

Recall that in ^U\, The Weyl asymptotic behavior of the eigenvalue counting function 
for the SQ case has been studied by Fukushima and Shima. Denote by p{x) the number 
of eigenvalues of A (taking the multiplicities into account) not exceeding x. According to 
their result, there exist positive constant c, C such that cx'^^l'^ < p(x) < Cx^^^^, for all x 
large enough, where ds = log9/log5 is the spectral dimension of SQ. In particular, 

< limmi p{x)x~'^^^'^ < limsup p{x)x~'^'^^'^ < oo. (1.1) 

Now what happens to the asymptotic behavior of the eigenvalue counting function on 
f2? A natural analogue of (11.11) holds which can be easily proved by first considering 
the asymptotic behavior of the eigenvalue counting function for each type of eigenvalues 
separately, then adding up them together. 

This is not the entire story. Recall that in the classical case. Suppose D is an arbitrary 
nonempty bounded open set in R" with boundary dD, then Weyl's classical asymptotic 
formula can be extended as follows: 

p{x) = (27r)-"c„|D|„x'^/2 ^ Q^^{n-l)/2^ 

as a; — > oo, where c„ depends only on n. See details in [261 ISHl [29]. The above remainder 
estimate constitutes an important step on the way to H. Weyl's conjecture [38] which 
states that if dD is sufficiently "smooth", then the asymptotic expansion of p{x) admits a 
second term, proportional to x^""^-*/^. Extending Weyl's conjecture to the fractal case, M. 
V. Berry [Sj H] conjectured that if D has a fractal boundary dD with Hausdorff dimension 
(which later was revised into Minkowski dimension in [6], [22]) dgn ^ {n — l,n\, then the 
order of the second term should be replaced hj dQci/2. See further discussion and a partial 
resolution of the conjectures of Weyl and Berry in M. L. Lapidus's work [22]. Hence it is 
natural to ask that is there an analogue result in SQ or Q setting. For SQ case, using a 
refinement of the Renewal Theorem, Kigami [18] showed that the remainder is bounded, 

p{x) = ^(logx)x'^s/2 + 0(l), 

where g is a. periodic function of period log 5. Note that this is consistent with the fact 
that the boundary of SQ consists of three points, hence has dimension zero. This was 
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refined by Strichartz in [36], where an exact formula was presented witli no remainder 
term at all, provided we restrict attention to almost every x. As for Q case. We will show 
that although we are unable to prove, it becomes apparent there is a second term of order 
log 2 / log 5 in the expansion of the eigenvalue counting function on Q from observing the 
experimental data. We note that our work deals with the case when the domain itself 
is fractal (and hence not open) also. The order of the second term should also has a 
close connection with the dimension of the boundary dQ duo to Weyl-Berry conjectures. 
Moreover, when consider a more general domain Q^, we will meet "drums with fractal 
membrane" with also fractal boundary. 

The paper is organized as follows. In Section 2 we briefly introduce some key notions 
from analysis on fractals and give a precise description of the Dirichlet and Neumann 
spectra of the Laplacian on SQ, which will be used in the rest of the paper. 

In Section 3, we first present the structure of the complete Dirichlet spectrum on 
Q before going into the technical details. We find an appropriate sequence of graph 
approximations for the fractal domain Q, and describe the structures of the corresponding 
discrete spectra of the successive difference operators on them. Accordingly, for each 
graph all the graph eigenvalues are also divided into three types, localized, primitive and 
miniaturized. By using an eigenspace dimensional counting argument, we show that they 
should make up the whole discrete spectrum. We also briefly describe how to relate 
the spectra of successive levels and how to pass the graph approximations to the limit by 
using spectral decimation for localized eigenvalues and weak spectral decimation for other 
types of eigenvalues. Then we list some conjectures concerning eigenvalue asymptotics 
(especially the existence of the second term of the expansion of the eigenvalue counting 
function), gaps in the ratios of eigenvalues and eigenvalue clusters which become apparent 
from observing the experimental data. At the end of this section, we present the structure 
of the Neumann spectrum on Q. 

In Section 4, we begin discussion of the discrete Dirichlet primitive eigenvalues. We 
will divide our discussion into symmetric case and skew-symmetric case. In each case, 
we will prove that for each level the primitive graph eigenvalues are the total roots of a 
high degree polynomial. And we describe the weak spectral decimation phenomenon by 
studying the relation between roots of consecutive polynomials. Moreover, we prove that 
the complete discrete spectrum is made up of the three types of eigenvalues as expected. 

In Section 5, we discuss the Dirichlet primitive eigenvalues on Q by passing the results 
of Section 4 on graph approximations to the limit. Since we can only use weak spectral 
decimation this time, some trivial results in SQ case become nontrivial and need to be 
proved in this section. 

In Section 6, we first prove that the whole Dirichlet spectrum on Q is made up of the 
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three types of eigenvalues following the basic idea of Fukushima and Shima's work. Then 
prove an analogue of Weyl's classical result on the eigenvalue asymptotics. The eigenvalue 
counting function p^{x) is shown to be of order x'^-^/^ as x — )■ oo where ds = log 9/ log 5 is 
the spectral dimension of SQ. Moreover, we also prove that the limit p^(x)/x^^/^ is not 
convergent. 

In Section 7, we give a brief discussion on how to deal with the Neumann spectrum. 
We will find a similar weak spectral decimation for primitive eigenvalues by establishing 
a relation between primitive symmetric (or skew-symmetric) graph eigenvalues with some 
high degree polynomials, but the proof is quite different from that in the Dirichlet case. 

We will also give a brief discussion on how to extend our method from Q to Qx with 
< X < 1 in Section 8. The purpose of this paper is to work out the details for one 
specific example. We hope this example will provide insights which will inspire future 
work on a more general theory. 



2 Spectral decimation on SQ 

First we collect some key facts from analysis on SQ that we need to state and prove 
our results. These come from Kigami's theory of analysis on fractals, and can be found 
in [151 [IS]. An elementary exposition can be found in [331 ESj. The fractal SQ will be 
realized as the limit of a sequence of graphs ro,ri,--- with vertices Vq C Vi C 
The initial graph Fq is just the complete graph on Vq = {qo,qi,q2}, the vertices of an 
equilateral triangle in the plane, which is considered the boundary of SQ. See Fig. 2.1. 
The entire fractal is the only 0-cell, which has Vq as its boundary. At stage m of the 
construction, all the cells of level m — 1 lie in triangles whose vertices make up Kn-i- 
Each cell of level m — l splits into three cells of level m, adding three new vertices to Vm- 




To Ti 

Fig. 2.1. The first 3 graphs, ro,ri,r2 in the approximation to the Sierpinski gasket. 
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We define the unrenormalized energy of a function u on by 

Em{u) = J2 (^(^) - ^(y))'- 

The energy renormalization factor is r = |, so the renormalized graph energy on Tm is 

and we can define the fractal energy £{u) = linim^oo £miu). We define domS as the space 
of continuous functions with finite energy. Then S extends by polarization to a bilinear 
form S{u, v) which serves as an inner product in this space. The energy S gives rise to a 
natural distance on SQ called the effective resistance metric on SQ, which is defined by 

d{x,y) = (min{£^(M) : u{x) = and u{y) = 1})^^ (2.1) 

for X, 1/ G SQ. It is known that d{x, y) is bounded above and below by constant multiples 
of \x — y|'°s(5/3)/i°g2^ where \x — y\ is the Euclidean distance. Furthermore, the definition 
(12. ip implies that functions on domS are Holder continuous of order | in the effective 
resistance metric. 

We let fj, denote the standard probability measure on SQ that assigns the measure 
3"™" to each cell of m level. The standard Laplacian may then be defined using the weak 
formulation: u G dorriA with Am = / if / is continuous, u G domS, and 



£{u, v) = — 




for all V G donioS, where donioS = {v E S : v\vo = 0} . There is also a pointwise formula 
(which is proven to be equivalent in [55]) which, for nonboundary points in = [J^ Vm 
(not in Vo) computes 

3 

Au{x) = - lim 5™A„m(x), 
where is a discrete Laplacian associated to the graph F^, defined by 

Amu{x) = ^ {u{y) - u{x)) 

for X not on the boundary. 

The Laplacian satisfies the scaling property 

A(uoFi) = \{Au)oFi 



and by iteration 

A(noF^) = i^(An)oF^ 

5™ 
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for F„ = F^, o o ■ ■ ■ o F^^. 



Although there is no satisfactory analogue of gradient, there is normal derivative 
dnu{qi) defined at boundary points by 



dnu{qi)= lim r - 



m— >oo 



the limit existing for all u G domA. The definition may be localized to boundary points 
of cells. For each point x G Vm\Vo, there are two cells containing x as a boundary point, 
hence two normal derivatives at x. For u G domA, the normal derivatives at x satisfy 
the matching condition that their sum is zero. The matching conditions allow us to glue 
together local solutions to Am = /. 

The above matching condition property follows easily from a local version of the 
following Gauss- Green formula, which is an extension of fl2.2p to the case when v doesn't 
vanish on the boundary: 



£{u, ~ ~ y {^'u)vd^ + vd^ 



Vo 

The local version of the Gauss-Green formula is 

£a{u,v) = — I {Au)vdii + ^ vdnU 



— I {Au)vdfi + vd, 



where A is any finite union of cells and £a{u, v) is the restriction of the energy bilinear 
form £{u, v) to A, which can also be defined directly by 

£a{u,v)= lim y^{u{x)-u{y)){v{x)-v{y)). 

m— >oo ' ' 

in A 

Now we come to a brief recap of the spectral decimation on SQ. Our goal is to find 
all solutions of the eigenvalue equation 

—Am = Am on SQ 

as limits of solutions of the discrete version 

-AmUm = XmUm OU Vm \ Vq. 

In the SQ case, we are lucky that we may take Um = u\v^, which is necessarily convenient 
for the spectral decimation. We should emphasize that this is not true for Q case. 

The method of spectral decimation on SQ was invented by Fukushima and Shima 
[To] to relate eigenfunctions and eigenvalues of the discrete Laplacian A^'s on the graph 
approximation F^'s for different values of m to each other and the eigenfunctions and 
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eigenvalues of the fractal Laplacian A on SQ. In essence, an eigenfunction on Tm with 
eigenvalue can be extended to an eigenfunction on r^+i with eigenvalue A^+i, where 
Am = /(Am+i) for an explicit function / defined by /(x) = x(5 — x), except for certain 
specified forbidden eigenvalues, and all eigenfunctions on SQ arise as limits of this process 
starting at some level mo which is called the generation of birth. This is true regardless 
of the boundary conditions, but if we specify Dirichlet or Neumann boundary condition 
we can describe explicitly all eigenspaces and their multiplicities. 
Denote the real valued inverse functions of /(x) by 4>±{x). That is 



(p±{x) = • (2.3) 

We describe the procedure briefly here. First, there is a local extension algorithm that 
shows how to uniquely extend an eigenfunction Um defined on Vm to a function defined 
on Kn+i such that the A-eigenvalue equations hold on all points of Kn+i \ K™- For SQ, 
the extension algorithm is: Suppose Um is an eigenfunction on with eigenvalue A^. 
Let Xm+i = 0±(Am)- Consider an m-cell with boundary points Xo,Xi,X2 and let yo,yi,y2 
denote the points in Kn+i \ Vm in that cell, with yi opposite Xj. Extend to a function 
Um+i on Vm+i by defining (for simplicity of notation, we drop the subscripts on u) 

Hy^) = 77^ r wT ^ , ^ = 0,l,2. (2.4) 

Then we have the following proposition taken from [35] . 

Proposition 2.1. Suppose Am+i 7^ 2, 5 or 6, and Am = /(Am+i)- Mm is a Am- 
eigenfunction of Am and is extended to a function u^+i on Kn+i by f l2.4p . then u^+i is 
a \m+i- eigenfunction of A^+i, Conversely, if Um+i is a Xm+i-eigenfunction of A^+i and 
is restricted to a function Um on Vm, then Um is a Xm- eigenfunction of Am- 

The forbidden eigenvalues {2, 5, 6} are singularities of the spectral decimation function 
/. It is "forbidden" to decimate to a forbidden eigenvalue. Because forbidden eigenval- 
ues have no predecessor, we speak of forbidden eigenvalues being "born" at a level of 
approximation m. 

Next we want to take the limit as m — )■ 00. We assume that we have an infinite 
sequence {Am}m>mo related by Am+i = 0±(Am) with all but a finite number of 0_'s. Then 
we may define 

A = ^ lim S'^A^. 

/ m^oo 

It is easy to see that the limit exists since 

0_(x) = -x + 0(x2) (2.5) 
5 

as X — !■ 0. Now suppose we start with a Amo-eigenfunction u of Am^ on Vmo, and extend u 
to K successively using (12. 4p . assuming that none of Am is a forbidden eigenvalue. Since 
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(12. 5 p implies Am = 0{-^) as m — )■ oo, it is easy to see that u is uniformly continuous on 
K and so extends to a continuous function on SQ. Moreover, it satisfies the A-eigenvalue 
equation for A. 

A proof in guarantees that this spectral decimation produces all possible eigen- 
values and eigenfunctions (up to linear combination). 

To describe the explicit Dirichlet and Neumann spectra, we have to describe all possible 
generations of birth and values for Am^, and describe the multiplicity of the eigenvalue 
by giving an explicit basis for the Amo-eigenspace of A^o- For each m, we have to add 
up the dimensions of eigenspaces with generation of birth mo < m, extended to Tm in all 
allowable ways. This total must be jlKn (Neumann) or jjV^ — 3 (Dirichlet), the dimension 
of the space on which the symmetric operator A^ acts. Now we give a brief description 
of the structure of the Dirichlet and Neumann spectra on SQ respectively. 

Dirichlet spectrum. 

We denote by V the Dirichlet spectrum of A on SQ and by the discrete Dirichlet 
spectrum of A^ on for m > 1. Due to the above discussion, we only need to make 
clear the spectrum Vm for each level m. There are two kinds of eigenvalues, initial and 
continued. The continued eigenvalues will be those that arise from eigenvalues of Vm-i 
by the spectral decimation. Those that remain, the initial eigenvalues, must be some of 
the fobidden eigenvalues by Proposition 2.1. 

In [30], it is proved that Vi consists of two eigenvalues 2 and 5 with multiplicities 1 
and 2 respectively, and for m > 2, the only possible initial eigenvalues in Vm are the two 
forbidden eigenvalues 5 and 6 with multiplicities ^"^ and ^"^"^ respectively. Hence we 
may classify eigenvalues into three series, which we call the 2-series, 5-series, and 6-series, 
depending on the value of Amo- The eigenvalues in the 2-series all have multiplicity 1, 
while the eigenvalues in the other series all exhibit higher multiplicity. Also, if A is an 
eigenvalue in the 5-series or 6-series, then S'^A is also an eigenvalue, corresponding to a 
generation of birth mo + m, with the same choice of 0± relations (suitably reindexed). 

Neumann spectrum. 

We impose a Neumann condition on the graph Tm by imagining that it is embedded 
in a larger graph by reflecting in each boundary vertex and imposing the Am-eigenvalue 
equation on the even extension of u. This just means that we impose the equation 

(4 - XmHq,) = 2u{Frq^+i) + 2u{F^q^-l) 

at Qi for 2 = 0,1, 2. Then the Neumann Am-eigenvalue equations consist of exactly jlKn 
equations in jj^n unknowns. Similar to the Dirichlet case, we also only need to make clear 
all the discrete spectra. The result is very similar to the Dirichlet spectrum, with only a 
few changes. We omit it. 
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It should be emphasized here that those eigenfunctions which are simultaneously 
Dirichlet and Neumann play an important role in the spectral analysis of SQ. Here 
we call them localized eigenfunctions since all of them have small supports. (Here the 
definition of localized eigenfunctions is slightly different from that of [21 [181 ESj for the 
convenience of further discussion for Q case.) Similar to V, to describe the structure of 
localized eigenfunctions, we only need to make clear the structure of all initial localized 
eigenvalues with generation of birth m for each value of m, which consists of 5-series and 
6-series eigenvalues. In fact, the multiplicity of a 5-series eigenvalue is pm(5) = 
with an eigenfunction associated to each m-level loop (a m-level circuit around an empty 
upside-down triangle in the graph r^). The eigenfunction u associated to each loop takes 
value on all m-level points not lying in that loop. Moreover, the support of u is exactly 
the union of all m-cells intersecting that loop. The multiplicity of a 6-series eigenvalue is 

(6) = ^ with an eigenfunction associated to each point x in Vm^i \ Vq. Each such 
eigenfunction u takes value on all points in Kn-i except x. Moreover, u is supported in 
the union of two (m — l)-level cells containing x. The existence of localized eigenfunctions 
is unprecedented in all of smooth mathematics. However, for a class of p.c.f. self-similar 
sets, including SQ, localized eigenfunctions dominate global eigenfunctions. See more 
details in [18]. 

3 The structures of Dirichlet and Neumann spectra 
on Q 

To give the reader an intuitive perception of the structure of the spectrum of A on f2 in 
advance, in this section we describe all Dirichlet or Neumann eigenvalues and eigenfunc- 
tions and their multiplicities on Q without proof. We will go to the technical details in 
the following sections. 

3.1 Dirichlet spectrum 

We begin with the Dirichlet case. Let S denote the Dirichlet spectrum of A on Q. We 
will consider three kinds of eigenfunctions, primitive, localized and miniaturized. 

The localized eigenfunctions are just a subspace of the localized eigenfunctions on SQ 
whose supports are disjoint from L (the line segment joining qi and ^2)- We let C denote 
the eigenvalues associated to them. These have generation of birth mo > 3 (the ones with 
mo = 2 all have supports intersecting L) and A^o = 5 or 6. 

Comparing to the SQ case, instead of the eigenfunctions associated to the 2-series 
eigenvalues, there is a type of global eigenfunctions, which we will call primitive eigen- 
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functions, which are sorted into symmetric and skew-symmetric parts according to the re- 
flection symmetry fixing go- We let and V~ denote the symmetric and skew-symmetric 
primitive eigenvalues associated to them respectively. And let V denote all of this kind of 
eigenvalues. They have multiplicity one. (An explanation of this will come later.) In fact, 
we call an eigenfunction u a symmetric primitive eigenfunction if it is symmetric under 
the refiection symmetry fixing go and also local symmetric in each cell F^{SQ) under the 
refiection symmetry fixing -F„,go with word w taking symbols only from {1,2}. For the 
skew-symmetric case, instead, we require u to be skew-symmetric under the reflection 
symmetry fixing go, but still local symmetric in small cells. Fig. 3.1. gives a symbolic 
picture of the symmetries, indicated by dotted lines. The primitive eigenfunction u (either 
the symmetric or skew-symmetric case) is unique and determined by the values denoted 
by {bo, bi, b2, ■ ■ ■) of u on vertex points (go, -Figo, F^qo, ■ ■ ■ ) by using the eigenfunction ex- 
tension algorithm described in (12. 4p . Due to the Dirichlet boundary condition, we always 
have bo = and lim„^_^ b^ = 0. We call (go, -Figo, F^qq, ■ ■ ■ ) a skeleton of fl since it plays 
a critical role in the study of primitive eigenfunctions. Another equivalent definition of 
the primitive eigenfunctions will be presented later by considering first primitive graph 
eigenfunctions then passing the approximation to the limit. 



bo 




Fig. 3.1. The first 4 level symmetries and the skeleton of ^. 



There is another type of eigenfunctions that we call miniaturized eigenfunctions. For 
each A G V~ there is a family of eigenfunctions with eigenvalue and multiplicity 2^ 
for A; = 1, 2, 3, ■ ■ ■ . To get such an eigenfunction, just take the A-eigenfunction u, contract 
it k times, place it in any one of the 2^ bottom cells of level k, and take value elsewhere. 
See Fig 3.2. The reason we can do this is that on the boundary point go the matching 
condition of u holds automatically. Let Ai denote all the eigenvalues associated to them. 

In section 6 we will prove that all eigenfunctions of A on 1] fall into one of these three 
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k^2 



Fig. 3.2. The first 2 level miniaturized eigenfunctions. 

types, and there are no coincidences of eigenvalues among different types. That is 

S = CUV U M (disjoint union). 

For the purpose of studying the structure of the spectrum S of A, the first thing we 
should consider is to describe the spectra of A^'s on the associated graph approximations. 
Similar to the SQ case, the fractal domain Q can be realized as the limit of a sequence 
of graphs Qm- More precisely, Vm > 1, let be a subset of Vm with all vertices lying 
along L removed. Let Qm be the subgraph of T^ restricted to V^. Denote by dflm the 
boundary of the finite graph Qm- It is easy to find that \ dflm and dQm approximate 
to Q and dQ as m goes to infinity respectively. See Fig. 3.3. We denote by Sm. the 
discrete Dirichlet spectrum of on Q^n for fri > 2. On Q^, the Dirichlet Am-eigenvalue 
equations consist of exactly tl(V^ \ dflm) equations in tKKri \ dflm) unknowns. We start 
from m=2 since there is no Ai-eigenvalue equation. For simplicity, let = ^(V^ \ dflm)- 
It is easy to check that a2 = 5, = 24, and more generally, 

qm+l 1 

am = - 

noticing that a„, = a„_i + 3™ + 2™"! - 3 ■ 2™-\ where 3"^ = tl(K^ \ Kn-i), 2"*-^ is the 
number of points lying on the bottom boundary of flm~i, and 3 ■ 2*""^ is the number of 
points in Vm\ Vm-i lying on L or dflm- 

Due to different types of eigenvalues of A, we should consider the associated different 
types of graph eigenvalues of A^. We now describe how to define Cm, Vm and M.m 
respectively. In fact, by the spectral decimation recipe, each localized eigenfunction u of 
A whose generation of birth mo < m can be restricted to flm to get a graph eigenfunction 
Um of Am with the Dirichlet boundary condition on dflm holding automatically. We 
call all this kind of graph eigenfunctions m-level localized graph eigenfunctions and all 
the associated eigenvalues are denoted by Cm,- However, we can not imitate this process 
to get the m-level primitive graph eigenfunctions since the Dirichlet boundary condition 
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^2 



Fig. 3.3. The first 3 graphs, ili, Q2, ^3 in the approximation to il. with inside points and 
boundary points represented by dots and circles respectively. 



would be destroyed if we do the similar restriction. But we can define m-level primitive 
graph eigenfunctions on Qm directly using the following way. We call an eigenfunction 
Urn of Am a m-level symmetric primitive graph eigenfunction if it is symmetric under 
the reflection symmetry fixing go and also local symmetric in Fw{SQ) fl under the 
reflection symmetry fixing Fy,qQ with word w taking symbols only from {1,2}. Denote by 
all the eigenvalues associated to them. V~ can be defined in a similar way. Let Vm 
denote all of them. The primitive graph eigenfunction Um (either the symmetric or skew- 
symmetric case) is unique and determined by the values denoted by {bo,bi,b2, ■ ■ ■ ,bm) 
of Um. on vertex points (go, -^i^o, F^qo, ■ ■ ■ , F^qo) by using the eigenfunction extension 
algorithm described in (12. 4p . Due to the Dirichlet boundary condition, we always have 
bo = bm = 0. We call (go, -Figo, F^qQ, ■ ■ ■ , F^q^) a skeleton of fi^. It also plays a critical 
role in the study of primitive graph eigenfunctions. Miniaturized graph eigenfunctions on 
VLra can be defined in a similar way by using miniaturization of skew-symmetric primitive 
graph eigenfunctions whose level strictly less than m. Denote by Aim all the associated 
eigenvalues. 

It is not difficult to describe all the localized graph eigenvalues in since they are 
almost the same as the SQ case. There are two kinds of eigenvalues in Cm-, initial and 
continued. The initial eigenvalues are 5 and 6. For the 6-eigenfunctions of A,„i on 
comparing to the 6-eigenfunctions of on F^, the only difference is those eigenfunctions 
whose support intersecting the boundary dVtm should be removed. A similar analysis 
shows that they are indexed by points in V^^i \ dVL„i-i. Hence the multiplicity of 6 is 
Pm(6) = (^m-i = ^"2"^ ~ 2™. Similarly, the 5-eigenfunctions of A^ on Vtrn are indexed by 
m-level loops except those loops touching dVL„i. Hence the multiplicity Pm(5) = Pm(5) — 
(1 + 2 + 22 + ■ ■ ■ + 2'"-2) = 3'"-'+! _ 2™-i. So that is the story for initial eigenvalues. The 
continued eigenvalues will be those that arise from eigenvalues of Cm-i by the spectral 
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decimation. Note that every eigenvalue A^-i of A^-i bifurcates into two choices of 
Am of Am by (12. 3p . except Am,~i = 6, which just yields the single choice Am = 3 since 
the other is a forbidden eigenvalue 2. We know that Pm_i(6) of Cm~i correspond to 
eigenvalue 6 of Am-i, while the remaining tl-Cm-i — Pm-i(6) of them correspond to other 
eigenvalues, leading to a space of continued eigenfunctions of dimension 2 • (jj£m-i — 
P™-i(6)) + P^_i(6) = 2 ■ tt^m-i - + 2—1. If we add to this p^{6) = - 2"^ 

and Pm(5) = 2 "^^ ~ 2™~i, we should obtain tl£m- Hence we have 

om— 1 1 om 1 om — 1 _i_ i 

tJLm = 2 ■ fll-m-i ^ K 2 H 2 H 2 

3™ + 1 

= 2 ■ tl£m-i H ^ 2 . 

Combining this with jj£2 = 0, we can easily get 

om+l 1 

= (m - 2) ■ 2" - 26 ■ 2'"-^ for m > 2. 

As for primitive graph eigenvalues Vm, things become more complicated. We consider 
"Pm ^iid respectively. We will show in the next section the spectral decimation recipe 
for this type of eigenvalues can not be used directly. In fact there is even not an analytic 
relation between elements in P+ (or V~) and elements in Pm+i (^^ "^m+i)- ^ rough but 
intuitive explanation of why does this "bad" thing happen is that the Dirichlet boundary 
condition will be destroyed when we use the eigenfunction extension algorithm (12. 4p to 
extend a Am-eigenfunction Um from Qm to flm+i or restrict a Am+i-eigenfunction Um+i 
from flm+i to Qrn- Howcvcr, a weak but useful relation between (or V~) and Pm+i (or 
V~^i) will be found in the next section, which will take the place of spectral decimation 
in the further discussion. We will prove that: For each m > 2, "P^ consists of = 
2™ + 2™"^ — 2 distinct eigenvalues with multiplicity 1, between and 6 strictly, denoted 
by Am,i, Am,2, ■ ■ ■ , Am,r™ in increasing order. Moreover, r^+i = 2rm + 2 and 

< Am+1,1 < 0-(Am,l), 

0-(Am,fe-l) < Am+l,fc < 0-(Am,,fc), V2 < A; < 

4>~{^m,rra) < Am+l,r,„ + l < 0+(-^m,r™), 
0+(Am,2r™+2-fc) < Am+l,fc < 0+ (Am,2r,„+l-fc) , Vr^ + 2 < < 2rm, 
4'+{^m.,l) < Am+l,2rm+l < 5, 
5 < Am+l,2r-m+2 < 6. 

Similar property holds for with rm replaced by Sm = 2™ — 2. In order to study the 
relation between P+ (or P") and Pm+i (or Pm+i)' introduce the following notations. In 
symmetric case, let 0_(Am,i) denote the (m + l)-level eigenvalue between and </)_(Am,,i)- 
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Let (f)-{Xm,k) denote the (m + l)-level eigenvalue between (f)_{Xm,k-i) and 0_(Am,fc) for 
each 2 < k < r^. Let (t>+{Xm,k) denote the (m + l)-level eigenvalue between 0+(Am,fc) 
and (f)+{Xm,k-i) for each 2 < k < r^- Let 0+(Am,i) denote the (m + l)-level eigenvalue 
between </)+(Am,i) and 5. Call this kind of (m + l)-level eigenvalues continued eigenvalues. 
There are another two (m + l)-level eigenvalues: one is between 0_(Am,r™) and 0+(Am,r™)5 
the other is between 5 and 6. Call these two (m + l)-level eigenvalues initial eigenvalues 
with generation of birth m + l. For the 2 level, all r2 = 3 primitive symmetric eigenvalues 
-^2,1; -^2,2 and A2,3 are called initial eigenvalues with generation of birth 2. We define the 
similar notations for skew-symmetric case in an obvious way with replaced by Sm- 
From this point of view, the continued primitive eigenvalues in P^+i (or V~_^i) will be 
those arise from eigenvalues in (or V~) by a (f)± bifurcation similar (but never equal) to 
(j)± bifurcation. We call this phenomenon weak spectral decimation, which will be proved 
playing a critical role in the study of the structure of primitive eigenvalues on Q in stead 
of spectral decimation. We should emphasize here that (j)± is not a real function relation. 
(It is just a notation for simplicity.) See the following diagram for the relation between 
and P^+i- The skew-symmetric case is similar. 

y \ ^/ \ ■■■ y \ 

The structure of Aim depends on the structure of all P^T's with A; < m by the definition 
of M-m- In fact, it is easy to check that 

m— 1 m— 1 

^Mra = S^'^'tl^fc" = Yl 2'^-''{2!' - 2) = (m - 3) ■ 2"^ + 4 for m > 2. 

k=2 k=2 

It will be proved that different types of eigenfunctions of on Q^n are linearly 
independent in the next section. Moreover, it is easy to check ]^Cm, ^Vm and jjA^m add 
up to tt(l^^ \ dVlm) since 

Qm+l 1 

^CmHVn.HMm = ^ (,„_2).2'"-26-2'"-3+r^+s^+(m-3)-2'"+4 = a^. (3.1) 

Hence we have the complete spectrum Sm of A^. It consists of the only above mentioned 
three types of eigenvalues. In Table 3.1, we list the eigenspace dimensions of all different 
types of eigenvalues in Sm for level m = 2, 3, 4, 5. 



level 






tJ-Pm 




tt<Sm 


m 


^"'^2'"^ (m 2) -2™ 26-2™-3 


2m _,_ 2"i-2 _ 2 


2™ - 2 


(m - 3) ■ 2™ + 4 


3m + l_l rnn + l 
2 ^ 


2 





3 


2 





5 


3 


6 


8 


6 


4 


24 


4 


37 


18 


14 


20 


89 


5 


164 


38 


30 


68 


300 
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Table 3.1. Eigeiispace dimensions of different types of eigenvalues in 5m- 

Next we want to take the limit as m — >■ oo. For L case, we assume that we have an 
infinite sequence of localized graph eigenvalues {Am}m>mo related by relations, with 
all but a finite number of 0_'s. Then we define 

A = ^ lim 5™Am. 

By successively using the eigenfunction extension algorithm fl2.4p from a Amp-eigenfunction 
u of on fimo; one can extend m to a localized eigenfunction of A on 1] associated to 
A. This method generates all the localized eigenfunctions £ similar to the SQ case. For 
case, we also assume that we have an infinite sequence of type graph eigenvalues 
{Am}m>mo related by 0± relations, with all but a finite number of 0_'s. Then we define 

A = - lim S'^Am. 

2 m— ^-oo 

We will also show the existence of the limit A. As shown before, now we can not use 
the eigenfunction extension algorithm. However, we can still prove that A is a P"*" type 
eigenvalue of A on f2 by a nonconstructive method. Furthermore, all P"*" type eigenvalues 
come in this way. This will be done in Section 5 by using the weak spectral decimation. 
The V~ and cases are completely similar to the case. Hence we get the complete 
spectrum 5 of A on 1). 



3.2 Spectral asymptotics, ratio gaps and clusters 

In Tables 3.2, 3.3, 3.4 and 3.5 we present the eigenvalues and their multiplicities in S^. 
for level m = 2,3,4,5, where we use A^^^, A~^,, A^.fc to denote the fc'th P"*", £ 
type eigenvalues respectively, and use M.mi}ymi k) denote the miniaturized eigenvalue 
generated from A~, ^. 

The following conjectures list some interesting phenomena we observed from the tables. 
Conjecture 3.1. Let pm{x) denote the eigenvalue counting function of Sm, i-e., 
Pm{x) = tt{A„, G S.m : A„, < x}. Then p„(0L"~'^H5)) = 3'^ - 2^ for k<m. 
Remark. Here 3'^ — 2^ is the difference between ak and at-i- 

This conjecture suggests that the bottom 3'^ — 2*^ eigenvalues of the Dirichlet spectrum 
of Vt should be generated from the bottom 3*^ — 2^ eigenvalues in Sm and the largest of 
these eigenvalues should be lim„^oo |5"'0^ ''''(5) = c5^ for the appropriate choice of c. If 
we define the Dirichlet eigenvalue counting function 

p^(x) = tt{A G 5 : A < x}. 
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then we have p^{ch^) = 3*^—2'^. This suggests an asymptotic growth rate p^(x) ~ 3;iog3/iog5 
as X — > oo. In analogy with the Weyl asymptotic law in SQ case, noticing high multiplicity 
of 0'^™~'^)(5) (which is ^'^ — 2*^"^) in Sm, using similar arguments we can get that the 
ratio p^{x) /x^°^^^^°^^ is bounded above and bounded away from zero, but non-convergent 
as X — oo. (This Weyl asymptotic law can be still proved without using Conjecture 3.1, 
by first considering the asymptotic law of the eigenvalue counting function for each type 
of eigenvalues separately, then adding up them together. See Section 6.) Moreover, of 
course, 

pf^j^^,-) ^ ^log3/log5 _ ^log2/log5 

along the sequence x = c5^. Hence, in analogy with the SQ case, we hopefully believe the 
following more precise conjecture. 

Conjecture 3.2. There exist two periodic functions gi{t) and g2(t) of period log 5, 
which are bounded above, bounded away from zero, and necessarily discontinuous at the 
value logc, such that 

p^{x) = ^i(loga:)x'°S=^/'°S^ + ^2(logx)x^°s2/log5 ^ ^^^log2/log5)_ ^3 2) 

Here besides the leading term 5'i(log Weyl's formula, the asymptotic 

second term of the eigenvalue counting function appears. This is very analogous to the 
conjectures of Weyl and Berry. 

Conjecture 3.3. There exist gaps in the ratios of eigenvalues from the Dirichlet 
spectrum S of A. That is, we can find infinitely many pairs of consecutive eigenvalues X, 
X' with y > c for some constant c > 1. 

Remark. In fact, in the discrete spectrum Sm., one can observe that gap appears 
above each for k < m. Moreover, there are also smaller gaps below miniaturized 

eigenvalues. 

In [5] it was shown that on SQ there exist gaps in the ratios of eigenvalues. The exis- 
tence of gaps is an interesting phenomenon in itself, but it also has important applications 
to analysis on fractals. See details in [5], [18], [36]. Thus it is of great interest to know 
whether similar phenomenon exists for fractals other than SQ. In fact [39] shows that 
this is the case for Vicsek set. Also [8] investigates this question for a variant of the SQ 
type fractal. 

Conjecture 3.4. In the spectrum Sm, between consecutive 5 and 6 type localized 
eigenvalues, there is exactly one and one V~ type eigenvalue ( except the case that the 
two consecutive eigenvalues are 0-(5) and 0+(6) = 3, where there is nothing in between). 

Conjecture 3.5. In the spectrum Sm, the number of distinct eigenvalues between 5 — e 
and 5 goes to oo as m ^ oo for any e > 0. 
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Remark. This means in S there exist eigenvalue clusters, that is, arbitrarily many 
distinct eigenvalues in an arbitrarily small interval. 

We say the spectrum S exhibits spectral clustering. Clustering does not occur on the 
SQ case. Experimental evidence suggests that it does occur on the pentagasket [1] and 
on the Julia sets [9j. It is proved that in [7] it also does occur on Vicsek set. 



m = 2 


eigenvalue Xm 


3 cm \ 
2 '-f ^m 


multi 


type 


m = 2 


eigenvalue Am 


3 cm \ 


multi 


type 


1 


aJi=1. 064568 


39.92 


1 


V+ 


4 


Aj3=5.472834 


205.23 


1 


V+ 


2 


Aj" 1=3.381966 


126.82 


1 


■p- 


5 


A2^^2=5.618034 


210.68 


1 


V- 


3 


A^2=4.462598 


167.35 


1 


V+ 











Table 3.2. The 2-level eigenvalues in 52 in increasing order. 



m = 3 


eigenvalue Am 


3 cm \ 
2 'J '^m 


multi 


type 


m = 3 


eigenvalue Am 


3 [rm \ 
2 ^ 


multi 


type 


1 


A+i=0.187518 


35.16 


1 




11 


A3 4=3.902230 


731.67 


1 


•p- 


2 


AJi =0.558733 


104.76 


1 


V- 


12 


A+g=4.517231 


846.98 
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V+ 


3 


A;^2=0-805532 


151.04 


1 


V+ 


13 


A3 5=4.803115 


900.58 


1 


V- 


4 


A3 2=1.247636 


233.93 


1 


■p- 


14 


A+7=4.946726 


927.51 


1 


V+ 


5 


A+3=l. 329287 


249.24 


1 


V+ 


15 


' A3,i=5 


937.50 


1 


C 


6 


A3 3=3.059152 


573.59 


1 


V~ 


16 


A =5.424059 


1017.01 


1 


r+ 


7 


A+^=3.075910 


576.73 


1 


v+ 


17 


A3 g=5.429135 


1017.96 


1 


V- 


8,9 


M3(A2_i)=3. 381966 


634.12 


2 


M 


18,19 


A^3{A-2)=5.618034 


1053.38 


2 


M 


10 


A+5=3.713736 


696.33 


1 


V+ 


20-24 


A3, 2=6 


1125.00 


5 


c 



Table 3.3. The 3-level eigenvalues in ^3 in increasing order. 



m = 4 


eigenvalue Am 
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v+ 
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2 


M 


40 
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V+ 


8 
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605.11 
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V+ 


41 
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1 


r+ 


9 
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611.81 


1 


V- 


42 
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1 


V- 


10 


AJ4=0.843591 


790.87 


1 


T- 


43 
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1 


V- 


11 


A+5 =0.857718 


804.11 


1 


V+ 


44 
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4300.11 


1 


v+ 


12 


A+g=0.965805 


905.44 


1 


V+ 


45 
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4439.70 
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p+ 
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999.09 


1 


T- 


46 
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1 


V- 
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M4(A-2)=l-247636 
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2 


M 
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2 


M 


16 
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V+ 


49 
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V+ 
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V- 
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V- 


18 
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1 


v+ 


51 
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1 


v+ 


19 
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1 


c 
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6 


c 
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2812.50 


5 
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58 


A:^ig=5.423778 
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v+ 
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M 
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M 


28 
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1 


V+ 
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M4(A^2)=5-618034 
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4 


M 
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3170.59 


4 


M 
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A4,5=6 


5625.00 


24 


C 


33 


A4,3=3.618034 


3391.91 


1 


C 













Table 3.4. The 4-level eigenvalues in ^4 in increasing order. 
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Table 3.5. The 5-level eigenvalues in ^5 in increasing order. 
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3.3 Neumann spectrum 

Next we give a brief discussion of the Neumann spectrum of A. Similar to SQ case, 
we want to impose a Neumann condition on the graph by extending functions from 
by even reflection, and imposing the pointwise eigenvalue equation at the boundary 
points in dQm, which now have 4 neighbors. Then the Neumann Am-eigenvalue equations 
consist of exactly equations in jjV^ unknowns. It is even convenient to allow m = 1, 
in which case there are three equations associated to the boundary dQi and no others. 
In particular, on Qi we find eigenvalues Ai = corresponding to the constant function, 
and Ai = 6 corresponding to the two dimensional space of functions satisfying u{qo) + 
u{Fiqo) + u{F2qo) = which can be split into an one dimensional symmetric space and 
an one dimensional skew-symmetric space under the reflection symmetry fixing go- For 
simplicity, let bm = tlKn- easy to check that 6i = 3, 62 = 10, and more generally, 

om+l 1 om+1 I 1 

bm = arn + mm = 2^+' + 2^ + 1 = ^ - 2^. 

We denote 5^ the Neumann spectrum of A on 1] and the Neumann spectrum of 
Am on Qm respectively. still consists of three types of eigenvalues, localized, primitive 
and miniaturized, denoted by , and respectively. And correspondingly, 
consists of three types of eigenvalues, denoted by and M.^ respectively. Moreover, 

{V^) can also be split into symmetric part V'^''^ {V^'^) and skew-symmetric part 
P~'^('P~'^) in the same sense as the Dirichlet case. 

The structure of localized (graph) eigenvalues is very similar to the Dirichlet case, with 
only a few changes: The 6-series has multiplicity increasing by 1, namely the eigenfunction 
associated to qo, while the 5-series is unchanged. Hence p^'^(6) = + 1 = ^"^"^ — 2"* 

and p^'^(5) = p^{5) = ^^^^^ - 2'"-\ Vm > 1, where p^'^(6) and p^'^(5) denote the 
multiplicities of the m-level initial eigenvalues 6 and 5 respectively. A similar discussion 
shows that 

^CZ = 2- I^Z-i - P^m-iiG) + pTiG) + Pr (5). 

Hence we have 

t^m = 2 ■ tt>C,„_i + 2 , 

which yields that 

om+l _ 1 

tj£^ = 2"*+^ - (m - 1) ■ 2"^ for m > 1, 

since = 0. 

The structure of primitive (graph) eigenvalues V^iV^) is also similar to the Dirichlet 
case. We consider the symmetric and skew-symmetric case respectively. In symmetric 
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case, we will prove that: For each m > 1, "P^' consists of 2"^ distinct eigenvalues with 
multiplicity 1, between and 6 with 0, 6 included, denoted by Am,i = 0, Am,2, ■ ■ ■ , ^m,2^ = 
6 in increasing order. Moreover, 

0+(-^m,2™+i-fc+l) < ^m+l,k < <P+i^m,2"^+^~k), 72"^ + 1 < < 2"^"*"^ — 1. 

Similarly, there is also a weak spectral decimation which relates "P^'^ ^m+i by intro- 
ducing the following notations. Let (f)-^{Xm^k) denote the (m + l)-level eigenvalue between 
0_(Am,fc-i) and 4>-{\m,k) for each 2 < k < 2™. Let 4>+{Xm^k) denote the (m + l)-level 
eigenvalue between (j)j^{Xm,k) and 0+(Am,fc-i) for each 2 < k < 2™. Call this kind of 
(m + l)-level eigenvalues continued eigenvalues. Hence the continued primitive eigenval- 
ues in V^'^i will be those arise from eigenvalues of V^'^ \ {0} by a (f)± bifurcation similar 
(but never equal) to bifurcation. There are another two (m + l)-level eigenvalues: one 
is 0, called zero eigenvalue, the other is 6, called initial eigenvalue with generation of birth 
m + 1. See the following diagram of eigenvalues in . 

6 

/ \ _ 

0-(6) 0+(6) 6 

<^_(/._(6) 9i+0-(6) ?i-0+(6) <^+0+(6) 0.(6) (/.+ (6) 6 

For skew-symmetric case, the only difference is that there is no zero eigenvalue in V^'^ . 
V~'^ consists of 2™ — 1 distinct eigenvalues between and 6, including 6, where 6 is an 
initial eigenvalue with generation of birth m and the others are continued eigenvalues arise 
from previous level eigenvalues by a similar weak bifurcation. We have now the following 
decimation diagram of eigenvalues in V~'^ . 

6 

^ \ ^ 

As for miniaturized eigenvalues, the structure of M.^ depends on the structure of all 
P^'^'s with A; < m in a completely similar way as the Dirichlet case. In fact, it is easy to 
check that 

m—l m—l 

= Yl '^'^'^V^k'^ = Yl 2"'-'' {2^ - 1) = (m - 2) ■ 2"^ + 2 for m > L 

k=l k=l 



i 

r+''' : 

i 
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It is easy to check ti£^, and tl-^m ^dd up to tlKn) since 

om+l 1 

^CZ + iV^ + ^Ml = ^ 2'"+i-(m-l)-2'" + 2"^ + 2'"-l + (m-2)-2™+2 = 6^. 

Hence we have the complete Neumann spectrum of A^- Then a completely similar discus- 
sion leads to the Neumann spectrum of A on VL. In Table 3.6, we list the eigenspace 
dimensions of all different types of eigenvalues in 5^ for level m = 1,2, 3, 4, 5. 



level 












m 


a^Y'i 2™+i (m l)-2™' 


2™ 


2™ - 1 


(m - 2) • 2™ + 2 


3™ + -' + ! r,m 
2 ^ 


1 





2 


1 





3 


2 


1 


4 


3 


2 


10 


3 


8 


8 


7 


10 


33 


4 


41 


16 


15 


34 


106 


5 


172 


32 


31 


98 


333 



Table 3.6. Eigenspace dimensions of different types of eigenvalues in 5^ . 
All the unproved details except the conjectures will be proved in the following sections. 




In this section, we work with m-level graph approximation fi^, m = 2, 3, 4 ■ ■ ■ . Denote 
by Vm the totality of primitive eigenvalues of the discrete Laplacian A^- In the following 
we use /^"^ to denote the n'th iteration of /, n > 1. We let f^^\x) = x. li w = f^'^\x), 
w is called a successor of x of order n with respect to /, and x is called a predecessor of 
w of order n with respect to /. 
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We begin with "P^, the symmetric eigenvalues in Vm- Let u„i be a Am-eigenfunction 
of Am with Am e "P^. Denote by (60, ^i, ^2) ■ ■ ■ , ^m) the values of Um on the skeleton 
(go, F'lQo, F'lQo, ■ ■ ■ , of fim where 60 = &m = by the Dirichlet boundary condition. 

See Fig. 4.1. Write aJ""^ the successor of A^ of order {m — i) with 2 < i < m for simplicity. 
Assume that none of X^^^^^s is equal to 2 or 5. (Later we will show this assumption 
automatically holds for any A^ G V^-) The eigenfunction extension algorithm (12. 4p gives 
the value of Um on the four {i + l)-level neighbors of -Fi^o for each 1 < i < m — 1, shown 
in Fig. 4.2. Hence the Aj-™]-eigenvalue equation at the vertex F^qo gives 



h-i 




Fig. 4.2. Values of Um on neighbors of Flqo. 



(A 9h , (14-3Ag)b. + (6-Aa)6._, 

iA-X\^i)b, = 2h+, + (^_;^M^^^_^(m.)^ ' Vl<^<m-1, (4.1) 

which can be modified into 

l{Xt:l)h^i + s{X^l)b, + riX^hk^i = 0, VI < z < m - 1, (4.2) 

with l{x) = X — 6, s{x) = (2 — a;)(4 — x)(5 — x) — (14 — 3x) and r{x) = —2(2 — 
x)(5 — x). Still from the eigenfunction extension algorithm, Um is uniquely determined 
by (&i, ■ ■ ■ 5 bm-i)- Here (61, ^2, ■ ■ ■ , ^m-i) can be viewed as a non-zero vector solution 
of either of the above two systems of equations consisting of m — 1 equations in m — 1 
unknowns. Hence the determinants of them should both be equal to 0. For simplicity, we 
are interested in the second determinant although comparing to the first one, it brings 
the possibility that A-'"'' {2 < i < m) could be 2 or 5, which should be removed. 
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The determinant associated to system (14. 2 p is a tridiagonal determinant, 



/(A, 



m — 1 y 



s(A(r)) 



Hence A^ should be a solution of the following equation 



/(/("-3)(a;)) r(/(™-3)(x)) 



/(/(x)) .(/(x)) r(/(x)) 
/(x) s(x) 



0. (4.3) 



Conversely, if A^ is a root of the polynomial qm{x) and none of f^'KXmVs with < 
z < m — 2 is equal to 2 or 5, then Am, € V^. Hence we are particular interested in all 
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the root x's of the polynomial qm{x) excluding those satisfying / 
< i < m — 2. We list some useful facts about the polynomial qm{x) 
Proposition 4.1. Let m > 2, then 

(1) gm(0) > 0; 

(2) qm{5) > 0; 

(3) QmiQ) < 0; 

(4) qmi<P^--'\5))<0; 

(5) qmicP'^T-'\2))>0; 

(6) gm+2(0L™"'^(3)) < and g3(3) > 0. 
Proof. (1) We will prove a stronger result, 

g^+i(0) > 20g„,(0) > 



2 or 5 for some 



(4.4) 



for m > 2. This can be proved by induction. It is easy to check that ^2(0) = 26 > and 
^3(0) = 556 > 20^2(0) by a direct computation. If we assume gm(0) > 20gm-i(0) > 0, 
then the expansion along the first row of gm+i(0) yields that 

g^+i(0) = 26g„(0) - 6 • 20gm_i(0) > 26g„(0) - 6gm(0) = 20g„,(0) > 0. 

(2) It is easy to compute that g2(5) = 1 > and ^3(5) = 6 > 0. For m > 4, 
g,n(5) = qm-M - 20g„„2(0) > by using 

(3) It is easy to compute that ^2(6) = — 4 < 0, ^3(6) = —3392 < ^2(6) < and 



gm(6) = s(/(™"2)(6)) ■ g^_i(6) - r(/(™-2)(6)) . /(/(-3)(6)) 



■ qm~2 



(6) 
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for m > 4 by the expansion along the first row of qm{Q)- 

Consider a polynomial defined by gi{x) = — r{f{x))l{x), it is easy to check 

that gi{x) > 1 whenever x < —6. In fact, we can write gi{x) = (2 — /(x))(5 — /(x))(4 — 
/(x) + 2(x — 6)) — (14 — 3/(x)) by substituting the expressions for s{f{x)), r{f\x)) and l{x). 
Noticing that 4— /(x)+2(a;— 6) = x^ — 3x— 8 > 46 and f{x) < whenever x < —6, we have 
9i{x) > 46(2-/(x))(5-/(a;))-(14-3/(x)) = 46(/(x))2-319/(x)+446. Moreover, since 
/(x) < —66 whenever x < —6, we finally have gi{x) > 46(— 66)^ — 319 ■ (—66) + 446 > 1. 

Then we can prove gm(6) < qm-i{6) < by induction. Suppose gm-i(6) < qm-2{6) < 
0. (This is true for m = 4.) Write g^(6) = ag„_i(6) + 6g,n-2(6) with a = s(/(™"2)(5)) 
and b = -r(/("'-2)(g)) . /(/(™-3)(6)). Noticing that m > 4, we have /("^-3)(6) < -6 
and /('"-2)(g) < 0. Hence a + b = ^i(/(™"^H6)) > 1 and 6 < 0. So by the induction 
assumption, we have 

gm(6) < aqm-i{6) + bqm-i{6) = (a + 6)5^-1(6) < gm-i(6) < 0. 

Hence we always have q'm(6) < for m > 2. 

(4) For simplicity, denote = qm{4>^~^~^\^)) ■ By direct computation, we have 02 = 
— 4 < and 03 ~ —92.10 < 0. We will prove a stronger result, am+i < lOom < 0, 
Vm > 2. It holds for m = 2. In order to use the induction, we assume Om+i < 10a,n < 0. 
An expansion of am+2 along the last row yields that 

«^+2 = .(0L"^+'H5))«™+i -r(0L"^)(5))/(0L"^+^H5))«^. 
Since 2 - (/.L"'^(5) > 0, 5 - 0L™^(5) > and (/>L"'+^^(5) - 6 < 0, we have 

= s(0L"^+')(5))a„+i-^r(0L"H5))/(0L™+'^(5))-(lOa^) 

< .(0L"^+^)(5))a„+i - lr(0(!"H5))/(0L™+^)(5))a^+i 

= [s{<P'"'^'\5)) - ^r(0L"^)(5))/(0L™+^)(5))]a^+i. 

Consider a polynomial 

92(3^) = s(x) r(f(x))l(x) = 14 + 9x + — x^ -x'^ H — x^. 

It is easy to compute that 

g',{x) = 9 - — X + — x^ - -x^' + x^ > 9 - —{<l>^'\5)) - -{<l>^'\5)r ^ AM > 
5 5 5 5 5 

whenever < x < 0^''(5). Hence g2{x) is monotone increasing in the interval [0,0^'* (5)]. 

Since < 0L"'"^^^(5) < 0L^^(5), we have ^2(0-""^^^ (5)) > ^2(0) > 10. Hence a^+2 < 
g2{(j)^r^'\5))am+i < lOam+i < 0. 
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The proofs of (5) and (6) are similar to that of (4). □ 

Now we discuss the possibihty of the roots of Qmix) satisfying (x) = 2 or 5 for some 
< i < m — 2. The following well-known basic algebra lemma should be useful. 

Lemma 4.1. Let g, h be two polynomials whose coefficients all belong to Q (the field 
of rational numbers), i.e., g,h E Q[a;]. If g is irreducible in Q[a;] and g, h have a common 
root in R, then g divides h in Q[x], i.e., all real roots of g belong to those of h. 

Lemma 4.2. Let x be a predecessor of 2 of order i with < i < m — 3. Then 
Qmix) = 0. Let X be a predecessor of 2 of order m — 2. Then qm{x) ^ 0. 

Proof. Firstly, let m > 3 and x be a predecessor of 2 of order i with < i < m— 3. Then 
= 2 and f'^'^^^\x) = 6. Substituting them into (14. 3p . noticing s(2) = r(6) = —8, 
s(6) = l{2) = -4 and r(2) = /(6) = 0, we get 



l(fii+i)(x)) r(/(*+i)(a;)) 

r(/«(a;)) 



-4 - 
-4 
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Secondly, let x be a predecessor of 2 of order m — 2. Then = 2. If m = 2, 

then X = 2. It is easy to check that x = 2 is not a root of g2(x). If m > 3, suppose x is a 
root of Qmix), then using the basic algebraic lemma Lemma 4.1, all roots of /*^™~^''(x) — 2 
are roots of qm{x). Noticing that 0^~^^(2) is also a root of f^'^~'^\x) — 2, we have 



,(m-2) 



(2)) = 0. But 



qn.{<P^r-'\2)) 



3(2) r(2) 
/(042)) s{<P42)) r(042)) 



/(0L'"-'^(2)) s(</.L™~'^(2)) 
s(0„(2)) r(0„(2)) 



s(2) 



■ Qm-li 



l{<P^^-^\2)) s{4>^!^-'\2)) 
^(_— 2)^2)), 



since r(2) = 0. By using Proposition 4.1(5), we get q„ 



,(m-2) 



(2)) < which contradicts 



to q„ 



, (m-2) 



(2)) = 0. Hence g™(x) ^ 0. □ 
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Lemma 4.3. Let x be a predecessor of 5 of order i with < i < m — 2. Then 
Qmix) ^ 0. 

Proof. Let x be a predecessor of 5 of order i with < i < m — 2. Then /'•*•' (x) = 5. 
Hence if x is a root of qmi^), then using Lemma 4.1, all roots of /'•'^(x) — 5 are roots of 
qm{x). Noticing that 0^''(5) is also a root of /'■'^(x) — 5, we have gm(0-''(5)) = 0. But 
9m (0^'' (5)) = qm{^) > by Proposition 4.1(2). More generally for < i < m — 2, 



/(5) s(5) r(5) 

/(0_(5)) s(0.(5)) r(0_(5)) 



/(0^^(5)) .(0^^(5)) 



s(/("^-2-^n5)) r(/("^-2-^)(5)) 



/(5) 



s(5) 
/(045)) s(045)) r(</)45)) 



/(0^^(5)) .(0«(5)) 



since r(5) = 0. Thus 



s(045)) r(0„(5)) 



/(5) .(5) 



/(0«(5)) s(0^^(5)) 



= g„„,(5)-g,+i(0^)(5))<O 

by the 2'nd and 4'th statements in Proposition 4.1. Hence VO < z < m — 2, we have 
proved gm(0-''(5)) ^ which yields a contradiction to gm(</'-^(5)) = 0. So qm{x) 7^ 0. □ 

From Lemma 4.2 and Lemma 4.3, for m > 3, the total unwanted roots of qm{x) are 
those predecessors of 2 of order i with < i < m — 3. q2{x) does not have any unwanted 
root. Hence to exclude them out, we define 



qm[x) 



(x-2)(/(x)-2)---(/(™-3)(x)-2)^ 



for m > 3, 
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and 

P2{x) = q2{x) = s{x). 

Pm{x) is still a polynomial from Lemma 4.2, although it looks like a rational function. 
Now we can say if is a root of the polynomial Pm{x), then G V^. Note that the 
degree of the polynomial qm{x) is 3 + 3 ■ 2 + ■ ■ ■ + 3 ■ 2™"^ = 3(2™"^ — 1) and the number 
of all the unwanted roots oi qm{x) is 1 + 2 + ■ ■ ■ + 2"*"^ = 2™"^ — 1 for m > 3 and for 
m = 2. Hence it is easy to check that the degree of Pm{x) is r„i = 2™ + 2*""^ — 2. The 
following is a list of some useful facts about the polynomial Pm{x). 
Proposition 4.2. (1) (-1)>^(0) > 0, Vm > 2; 

(2) p2{5) > and (-1)™-Vm(5) > 0, Vm > 3; 

(3) p2{6) < and (-1)>„(6) > 0, Vm > 3. 

Proof. It can be checked by a direct computation when m = 2. When m > 3, noticing 
that by the definition oi pm{x), 

P™(0) = 7^, P„(5)- 



and 

Pm(6) 



gm(6) 



(6-2)(/(6)-2)---(/(-3)(6)-2)- 

Using Proposition 4.1(l)-(3), we get the desired result. □ 

We now present a more precise result about the distribution of the roots oipm{x) and 
show an useful relation between roots of two consecutive polynomials. 

Lemma 4.4. For each m > 2, Pm{x) has distinct real roots satisfying 

< \m,l < Am,2 < ■ ■ ■ < Am,r-m-l < 5 < A^.r^ < 6. 

Moreover, (-l)'"+'=-Vm+i(</'-(A™,fc)) > and (-l)'"+*^p„+i((^+(A„,fc)) > 0, VI < < r^. 
Proof. We prove it by using the induction on m. 

When m = 2, P2{x) = s{x) has 3 distinct roots: A2,i ~ 1.0646, A2,2 ~ 4.4626 and 
A2,3 ~ 5.4728 by a direct computation. 

Let A be one of A2,fc's, then p2(A) = 0, i.e., s(A) = 0, and p3(0_(A)) = '^^^'^^|^2 ~ 
^i'l>'W-<in^-^)i5-^) by using s(A) = 0. Since < A < 6, we have 0„(A) - 2 < and 
0_(A) — 6<0. Hence p3(0_(A)) ~ (2 — A) (5 — A) where "~" means both sides of "~" have 
the same signs. Similarly, p3(0+( A)) = ^('^+('^)^6K2-a)(5-a) ^3(0,^ (A)) ~ — (2 — A)(5— A). 

Hence < A2,i < 2 yields that P3(0_(A2,i)) > and p3(0+(A2,i)) < 0; 2 < A2,2 < 5 
yields that P3(0-(A2,2)) < and P3(0+(A2,2)) > 0; A2,i > 5 yields that P3(0-(A2,3)) > 
and p3(0+(A2,3)) < 0. So our lemma holds for m = 2. 

We now assume our lemma holds for m, and prove it for m + 1. 
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Noticing that from Proposition 4.2, we have pm+i(0) ~ (—I)'" ^, Pm+i(5) ~ (— 1)"^ 
and pm+i{Q) ~ (—1)'""-'^. Hence if we write 

0, 0_(Am,i), 0_(A„^,2), ■ ■ ■ , 0-(Am,r„), 0+(Am,r„), ■ ■ ■ , 0+(Am,2), (f)+{\m,l), 5, 6 (4.5) 

in increasing order, then the values of Pm+i on them have alternating signs by the induction 
assumption. Hence there exist at least 2rm + 2 = r^+i distinct roots of pm+i{x), with 
each located strictly between each two consecutive points in (14. 5p . Moreover, these are 
the totality of the roots of pm+i{x) since the degree of pm+i{x) is also Vm+i- Hence we 
can write them in increasing order: 

< Am+1,1 < Am+1,2 < ■ ■ ■ < Am+l,r„i+i-l < 5 < A^+i^^^^j < 6. 

Now we study the signs of Pm+2(0±(Am.+i,fc))'s- Let A be one of Am+i,fc's, thenpm+i(A) = 
0. Moreover, 

/, gm+2(0-(A)) 
Pm+2(0-(A)) = - 



(04A)-2)(A-2)---(/(-2)(A)-2) 
s(0_(A))g^+i(A) + 2(0_(A) - 6)(2 - A)(5 - A)g^(/(A)) 
(0„(A)-2)(A-2)---(/(™-2)(A)_2) 

by using the expansion of qm+2{4>-W) along the last row. Since Pm+i(A) = 0, we have 
?m+i(A) = 0. Hence 

. , ^ 2(0-(A)-6)(2-A)(5-A)g^(/(A)) ^ -2(04A) - 6)(5 - A)p^(/(A)) 
p„,+2{<P-{^)) ^^_^A) _ 2)(A_2) ■■■(/(— 2)(A)- 2) 0-(A)-2 

Since < A < 6, we have 0_(A) — 2 < and 4>-{X) — 6 < 0, hence 

p™+2(0-(A))~(A-5)p^(/(A)). 



Similarly, 



and 



„ . . _ -2(0+(A) - 6) (5 - A)p^(/(A)) 

P^^2{M>^)) 0^(A)-2 



p„+2(0+(A))~(5-AK.(/(A)). 



When A = A^+i,!, we have < A < 0„(Am,,i), hence < /(A) < A^,!- Noticing 
that Am,i is the least root of Pm{x) and > by the induction assumption, we have 
PmifW) ~Pm(0) ~ (-1)". Hence p^+2(0_(A„+i,i)) ~ (-1)™+^ andp„+2(</'+(A„+i,i)) ~ 
(-1)™ since Xm+i,i < 5. 

When A = \m+i,k with 2 < k < r^, we have 0_(Am,fc-i) < A < (p^{\m.,k), hence 
Am,fc-i < /(A) < Xm,k- Noticing that Pm(Am,,fc-i) = and Pm.(O) ~ (-1)™, we have 
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Pm(/(A)) ~ by using the induction assumption. Hence 

(-1)™+^ and p™+2(0+(A„+i,fc)) ~ (-1)™+^-+! since A^+i.^ < 5. 

When A = \m+i;rm+i^ we have 0_(Am,r.„) < A < 0+(Am,,.„), hence /(A) > Am,r,„. 
Noticing that \m,rm is the last root of Pmi.^) and Pm(0) ~ (—1)"*, we have PmifW) ~ 
^_l^m+r,n by usiug the iuductiou assumption. Hence pm+2(0-(Am+i,r„+i)) ~ 
and Pm+2(0+(Am+i,^^+i)) ~ (-1)"'+^- since A^+i^^^+i < 5. 

When A = A^+i^^ with r,„+2 < A; < 2rm, we have (l)+{\ra,rm+i-k) < A < (j)+{\ra,r^+^-k-i), 
hence Am,^^_^^_fc_i < /(A) < A„,,.^_^^_fc. Noticing that Pm(A„,r^^j_fc_i) = and Pm{0) ~ 
(-1)"", we have Pm(/(A)) ~ ^ (_i)m+fc-i ^y using the induction as- 

sumption. Hence p„+2(0-(A™+i,fc)) ~ (-1)™+*^ and p„+2(0+(Am+i,fc)) ~ since 

Am+l,fc < 5. 

When A = \m+i,2rm+iy we have (f>^{Xm,i) < A < 5, hence /(A) < A^^i. So we have 
Pm(/(A)) ~Pm(0) ~ (-l)"". Hence pm+2(0-(A„+i,2r™+i)) ~ (-1)™+^ andpm+2(0+(A„+i,2r™+i)) - 
(-1)™ since A„+i,2r„,+i < 5. 

When A = Am+i,2r„+2, we have 5 < A < 6, hence /(A) < 0. So we have Pm,(/(A)) ~ 
Pm{0) ~ (-l)"". But now A > 5, hencep„+2(0-(A„,+i,2r„+2)) ~ (-1)™ andpm+2(0+(Am+i,2r„+2)) 

Hence we have proved (-l)"+^+^~Vm+2(0-(A™+i,fc)) > and {-l)"'+^+''p^+2{<|)+{K^+l,k)) > 
0, VI < /c < Vm+i- So our lemma holds for m + 1. □ 

Thus by Lemma 4.4, in particular the proof of Lemma 4.4 and the fact that each root 
of Pm{x) belongs to V^, we have the following result: 

Lemma 4.5. For eachm > 2, P+ consists of at least distinct eigenvalues satisfying 



< Xm,l < Ki,2 < ■ ■■ < Am,r™-1 < 5 < X„^,rm < 6- (4.6) 



Moreover, 



< A,„+i^i < 0_(Am,i), 

0_(Am,fc_i) < Xm+i,k < 0-(Am,fc), V2 < A; < r^, 

0-(Am,r„J < Am+l,r„+l < 0+(Am,r™), 

0+(Am,2r™+2-fc) < A„,+i,fc < 0+ ( Am,2r„+i-fc) , Vr„^ + 2 < < 2rm, (4.7) 

0+(Am,l) < Am+i,2r„+l < 5, 

5 < Am+l,2r„+2 < 6. 



Remark. The third inequality in (14. 7p can be refined into 2 < Am+i,r„i+i < 0+(Am,rm)- 
5*66 details in Theorem A in Appendix. 
Moreover, we have 
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Lemma 4.6. Let Am be a root ofpmix), Um a primitive Xm- eigenf unction on Qm, o,nd 
{bo,bi, - ■ ■ , b^) the values of Um on the skeleton of Qm- Then bi ^ and ^ 0. 

Proof. Without loss of generality, assume m > 3. We still use A^^'"'' to denote the 
successor of Am of order (m — i) with 2 < i < m. From the definition of pm{x), A^-"^^ 2 
or 5, for each 2 < i < m. From the discussion in the beginning of this section, the vector 
[bi, 62? ■ ■ ■ ) bm-i) can be viewed as a non-zero vector solution of system (14. 2p of equations. 

Suppose bm-i = 0. Then (61, 62, ■ " " ? &m-2) can be viewed as a non-zero vector solution 
of the system of equations consisting of the first (m — 2) equations of (14 ■2p in (m — 2) 
unknowns. Hence the determinant of this system gm-i(A^i) should be equal to 0. Thus 
X^m-i is a root of pm-i{x) since its all successors ■ ■ ■ , X^^l^ do not take value from 
{2,5} obviously. Then Lemma 4.4 says that neither of 0±(A^]^) should be a root of 
Pm{x). This contradicts to Pmi^m) = since A^ is equal to either of 0±(A|,'^]^). Hence 
bm~i 7^ 0. 

On the other hand, if bi = 0, then by substituting it into (14. 2p . noticing that none of 
Xl"^^'s is equal to 2 or 5, we can get 62 = 0, ■ ■ ■ , fem-i = successively, which contradicts 
to bm-i ^ 0. Hence 61 7^ 0. □ 

Next we give a brief discussion of the skew-symmetric case. It is very similar to 
the symmetric case. Let be a Am-eigenfunction of with A^ G V~. Denote by 
(607 ^1, ^2, ■ ■ ■ ,bm) the values of Um on the skeleton of Qm where bo = bm = by the 
Dirichlet boundary condition. Write A^-™'' the successor of A^ of order (m — i) with 
2 < i < m. Comparing to the symmetric case, the eigenvalue equations at the vertex 
F^go's are unchanged except the one at FiQq, since now the values of Um on the four 
2-level neighbors of FiQq are modified as shown in Fig. 4.3. Hence we still have the same 



u 




Fig. 4.3. Values of Um. on neighbors of Figo- 
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modified eigenvalue equation 

l{Xt:l)b^^l + + r(AH)6,+i = 0, \/2 < ^ < m 

while the first equation in f l4.2p is replaced by 



with s'(x) = (4 — a;)(5 — a;) — 1 and r (x) = — 2(5 — x). Now we assume A2'"'' 7^ 5 and none of 
A^^^'s is equal to 2 or 5 for 3 < i < m. (Later we will show this assumption automatically 
holds for any A^ G V^.) Then by the eigenfunction extension algorithm, is unique 
and determined by its values on the skeleton of Qm- Using similar discussion, A^ should 
be a solution of the following equation 

/(/(— s(/(™-3)(a;)) r(/(™"3)(x)) 



qm[x) 



l{f{x)) s{f{x)) r{f{x)) 
l{x) s{x) 



(4.8) 



instead of qm{x) = in the symmetric case. Hence if Am is a root of qm{x), f^"^~^\\m) 7^ 5 
and none of (Am)'s is equal to 2 or 5 for < i < m — 3, then A^, G V^. Similar to 
Proposition 4.1, we have 

Proposition 4.3. (1) qm{0) > 0, Vm > 2; 

(2) 52(5) < and q^{5) > 0, Vm > 3; 

(3) 52(6) > and q„^{6) < 0, Vm > 3. 

Proof. The first two statements follow from a very similar argument in the proof of 
Proposition 4.1. We only need to prove the third one. 

It is easy to check that $2(6) = 1 > and $3(6) = —436 < by a direct computation. 
For m > 4, an expansion along the first row yields that 

qm{G) = s{f^"'-'\6))qmMG) + 2(5 - f^^"'\6)){f^^-'\6) - 6)5^-2(6). 

Recall that in the proof of Proposition 4.1(3), we have proved that gm_i(6) < ^^-2(6) < 0. 
Hence 

qUQ) < mf^'^~'\Q)) + 2(5 - /(™"2)(6))(j(™-3)^Q) _ 6))g^„i(6), 

noticing that /(™-2)(6) < f^^^-^^G) < -6. An easy calculus shows that s(/('"~2)(g)) + 
2(5 - /(— 2)(6))(/(— 3)(g) _ 6) > 1, hence ^^(6) < g™_i(6) < 0. □ 

Similar to the symmetric case, the following two lemmas focus on the possibility of 
the roots of qm{x) satisfying /^''""^^(x) = 5, or /^*^(x) = 2 or 5 for some < z < m — 3. 
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Lemma 4.7. Let x be a predecessor of 2 of order i with < i < m — 3. Then 
Qmix) = 0. 

Proof. If < i < m — 3, the proof is the same as that of Lemma 4.2. So we only need 



to check the i = m — 3 case. In this case, / 



(m-3) I 



X) 



2 and = 6. Substituting 



them into (gS]), noticing s(2) = -8, l{2) = -4, r(2) = 0, s(6) = 1 and f(6) = 2, we get 



qm[x) 



r(/('"-3)(a;)) 



1 

-4 



-8 



0. □ 



Lemma 4.8. Let x be a predecessor of 5 of order i with < i < m — 2. Then 
qm{x) ^ 0. 

Proof. Similar to the proof of Lemma 4.3, we only need to prove qm{4'^^{^)) 7^ 0. A 
similar argument yields that gm(0^H5)) = ^^(5) and gm(0^^(5)) = gm-i(5) ■ gi+i(0^^(5)) 
for < i < m — 2. Combined with Proposition 4.1(4) and Proposition 4.3(2), it follows 
the desired result. □ 

Hence the total unwanted roots of qm{x) consist of those predecessors of 2 of order i 
with < 2 < m — 3 for m > 3 and q2{x) does not have any unwanted root. This is exactly 
the same as the symmetric case. To exclude them out, we define 



Pm{x) 



qm[x) 



{x-2){f{x)-2)---{f(^-^){x)-2y 



for m > 3, 



and 



P2ix) = q2[x) = s(x). 

These polynomials play a very similar role to Pm{x)^s in the symmetric case. It is easy to 
check that the degree of Pmix) is Sm = 2"^ — 2, since the degree of the polynomial qm{x) 
is3 + 3- 2 + -- - + 3- 2™-3 + 2 ■ 2'"^2 ^ 3(2"^-2 - 1) + 2™-\ and the number of all the 
unwanted roots of ^^(2;) is 1 + 2 + ■ ■ ■ + 2™"'^ = 2"^"^ — 1 for m > 3 and for m = 2. 
The following is a list of some facts about Pm{x) similar to Proposition 4.2, which can be 
easily get from Proposition 4.3. 

Proposition 4.4. Let m > 2, then 

(1) (-1)™F^(0) > 0; 

(2) (-l)'"-%(5) > 0; 

(3) {-lTpm{<o) > 0. 

Then following a similar argument, the results in Lemma 4.5 and Lemma 4.6 still hold 
with Pm{x) and replaced by V~, Pm{x) and respectively. 

Hence we have found distinct eigenvalues in P+ and Sm distinct eigenvalues in V~. 
We will show these eigenvalues are the totality of Vm- To prove this, the following lemma 
is needed. 
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Lemma 4.9. LetV^'* andV^'* be the sets of total roots of pm{x) and pm{x) respec- 
tively. Let Ai^ be the set of miniaturized eigenvalues generated by V^"* with 2 < k < m. 
Let Cm denote the set of m-level localized eigenvalues. Then all eigenfunctions associated 
to these eigenvalues are linearly independent. 

Proof. Without loss of generality, assume m > 3. It is easy to check that for each m- 
level localized eigenfunction m^, it must be on dVLm-i- Lemma 4.6 says that each m-level 
primitive symmetric Am-eigenfunction m^'^ with G V^* must be a non-zero constant 
on dVtjn-i \ {q'o} and be a non-zero constant on dVLi \ {q^}. The skew-symmetric analog of 
Lemma 4.6 says that each m-level primitive skew-symmetric Am-eigenfunction with 
Am G "P"'* must be a non-zero constant on each symmetric part of dVLm-i \ {go} under 
the symmetry fixing go? and take non-zero value on FiQq and F2qQ only different in signs. 
From the construction of the miniaturized eigenfunctions, for each m-level miniaturized 
eigenfunction u-^ with eigenvalue in M.^, must take non-zero value on a subset of 
dVLm-i \ {go} and be on dVti. These observations implies the linearly independence of 
eigenfunctions among different types. □ 

Hence we have 

Lemma 4.10. For each m>2, = P+ and P^-* = P^. 

Proof. It follows directly from Lemma 4.5 and its skew-symmetric analog. Lemma 4.9 
and the eigenspace dimensional counting formula (13. ip . □ 

By this lemma, it is easy to see that the assumptions we made before on symmetric 
and skew-symmetric m-level primitive eigenvalues hold automatically. 

Next we will prove each primitive eigenvalue A G Vm has multiplicity 1. For this 
purpose, we need the following lemma. 

Lemma 4.11. For each m > 2, 7^+ fl P+.^.i = 0. 

Proof. For m = 2, it can be checked by a direct computation. In order to use the 
induction, we assume that H V^+i = and will prove V^+i H = ^■ 

Suppose there is a A G V^+i H Vm+2- Then Pm+i(A) = Pm.+2(A) = (hence gm+i(A) = 
gm+2(A) = 0). Moreover, none of /^*^(A) (0 < z < m) is equal to 2 or 5. 

The expansion along the first row of gm+2(A) gives 

g™+2(A) = s(/(-)(A))g™+i(A) - r(/(™)(A))/(/(™-i)(A))gm(A). 

Noticing gm+i(A) = gm+2(A) = 0, we have 

r(/(-)(A))/(/('"-^)(A))g„,(A) = -2(2 - /("^)(A))(5 - f'''^\mf^"'-'\\) - 6)g„(A) = 0. 

Hence qm{\) = or /(™-i)(A) = 6, since /("^^(A) ^ 2 or 5. 

If gm(A) = 0, then A G V^, hence H V^^i 7^ 0. This contradicts to our induction 
assumption. 
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Hence we have /("^~^^(A) = 6, i.e., /'■"^"^^(A) = 3. Noticing that A is also a root of 
Pm+i{x), Lemma 4.1 says that 0^~^^(3) is a root of Hence Pm.+i(0^~^''(3)) = 0, 

which contradict to Proposition 4.1(6). Hence such A can not exist. So we get the desired 
result. □ 

Then we can prove: 

Lemma 4.12. For each m > 2, P+ fl V~ = 0. 

Proof. For m = 2 or 3, it can be checked by a direct computation. Let m > 4. Suppose 
there is an eigenvalue Am € fl V~. Then by Lemma 4.10, Pm{^m) = Pmi^m) = 0. For 
each 2 < i < m, denote by A^-'"'* the successor of Am of order {m — i). Obviously we have 
Qmi^m) = Qmi^m) = and Aj-™'' 7^ 2 or 5 for 2 < i < m. Furthermore, by Lemma 4.11, 
we have Pm,_i(Am) 7^ 0, hence gm-i(Am) 7^ 0. 

Using the expansions of gm.(Am) and gm(Am) along their first rows respectively, we have 

s(Af ))gm-i(Am) -r(Af ))/(Af ))gm-2(Am) = 

and 

?(Af ))gm_i(Am) -?'(Af))/(Af ))gm_2(Am) = 0. 

Hence, the vector {qm-i{^m), (lm-2{^m)) can be viewed as a non-zero solution of the system 
of linear equations, 

.(Ar))x-r(Ar))/(Af))y = 
^ ?(Af))x-f(Ar))/(Af))y = 0. 

Thus 

s{X^r^) 2(2-Af))(5-Af))(Af)-6) 
s{X^r^) 2(5-Af^)(Af)-6) 

Since A^™^ ^ 5, we have A^""^ = 6 or ^(A^'"^) = (2 - A^'"^)S'(A^™^). By Substituting the 
expressions for s{x) and s'(x), we get Ag"^'* = 6 or Ag™"* = 6. Hence we have Ag™"* = 3 or 
Af) = 3, i.e., /('"-^HArn) = 3, or /(™-^)(Am) = 3. 

Noticing that Am, is a root of qm{x), by using Lemma 4.1, we can see that either 
</)L"'"'^(3) or 0L™-^^(3) is a root of gm(a;), i.e., gm(0-""'H3)) = or gm(0-""^^(3)) = 0. An 
expansion of qm{4>^-^~^\'^)) along the first row yields that 



0. 



,(0(!"-')(3)) = .(/(^)(3))gm-i(0L'"-'^(3)) = 848gm-i(0L'"-'^(3)) 



since /(/(3)) = 0. Hence we have either qm{4>^-^ ^\'^)) = or gm-i(0^ ^'*(3)) = 0. By 
Proposition 4.1(6), this is impossible. Hence such Am can not exist. So ^'Pm = 0- n 
We summarize what we have proved: 

Theorem 4.1. For each m > 2 , consists of Vm distinct eigenvalues satisfying 

< Am,,l < Am,2 < • ■ ■ < Am,rm-1 < 5 < Am,r™. < 6. 
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A relation between and V^j^i is shown in i \A.7\i . Similar properties hold for with 
Tm replaced by Sm- Each G Vm has multiplicity 1. Moreover, the spectrum Sm of A^^ 
on Qrri satisfies 

where the union is disjoint. 



5 Primitive Dirichlet eigenvalues of A 

Having found the primitive Dirichlet eigenvalues and eigenfunctions for A^, it is natural 
to believe that the primitive Dirichlet eigenvalues of A can be obtained in the limit as m 
goes to infinity. This is true for the spectrum in SQ case, benefiting from the spectral 
decimation method and the eigenfunction extension algorithm. Our goal in this section is 
to extend this recipe to Q case by instead using the weak spectral decimation introduced 
in Section 3. Comparing to the SQ case, our method is more based on estimates. We will 
focus on the symmetric case, since the skew-symmetric case can be got by using a similar 
discussion. We use the 0± notations introduced in Section 3. Recall that if am, Pm are 
two consecutive eigenvalues in "P^ with ct^ < (3m, then we always have 

</)_(am) < 4>-{f3m) < (p-{(3m) and </)+(/3„) < (j)+{f3m) < (l>+{(^m), (5.1) 

and if Pm is the least eigenvalue in then instead we have 

< < (f)-{/3m) and (f)+{Pm) < < 5. 

Let niQ > 2, Xm^ be a mo-level primitive symmetric eigenvalue, {Xm}m>mo be an 
infinite sequence related by A^+i = 0-(Am) or 0+(Am), Vm > mo, assuming that there 
are only a finite number of 0+ relations. Call the minimum value mi, such that Vm > mi, 
Am+i = 4>-{Xm), the generation of fixation of the sequence {Xm}m>mo- In all that follows 
in this section, we always use {Xm}m>m.o as such a sequence without specifical declaration. 

The first fact about this sequence is: 

Lemma 5.1. lim-m^oo 5™Am exists. 

Proof. Without loss of generality, assume A^i < 5, otherwise, we could choose mi = 
mi + 1 and use m-i to replace mi in the following proof. 

Let m > mi then '^'"+^ = '^-^^"^'> < '?^-(^'") = <P-i^m.) _ i Since 

— ' Am Am — A„i (/)_(A,„)(5-</)_(Am)) 5 — 0_(Am)' 

< A„ < 5, we have < 0_(A^) < 2, hence < |. Thus Y.m>mi ^rn < oo. 

Fiirthprmnrp ^'""^^^'^+1 — c; Am+i ^ 5 _ i i (/'-(Am) NntiHncrthat V 
rurxnermore, ^^^^^ ~ ^ x„, - 5-(/._(Am) ~ -^^s-^-iAm) • ^^Q^^'^^^S ^^^^ Z^m>mi 5-0_( 

1 Y.m>rm (P-i^m) < | Em>mi < OO sluce < 1 wheucver < X < 5, wc get that 
nm>mi converges. Hence lim„_^oo S'^A™ exists. □ 
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The following is an estimate of the difference between 0_(Am) and 0_(Am) for Am in 
the sequence {Xm}m>mo- 
Proposition 5.1. 

J] 5"^(^_(Am)-</)_(Am)) <oo. 

m>mi 

In particular, limm-i-oo 5™'(0_(Am) — 0-(Am,)) = 0. 

Proof. Without loss of generality, assume A^i < 5. From Lemma 5.1, we have 
Em>™,(5'"+'Am+i - S'^Am) < oo. Hence 

m+l ^mj 

m>mi 

J2 5™(50_(Am) - 0_(Am)(5 - 0_(Am))) 

m>mi 

= E (5""''(0-(A-) - 0-(A-)) + 5™(0-(Am))') < oo. (5.2) 

m>mi 

Since < < 1 whenever < x < 5, we have 5™'(0„(Am))^ < S'^A^^. Still 

from Lemma 5.1, we have Am = 0{-^), hence 5™(0_(Am))^ < for some constant c. 
Thus Em>mi 5'"(0_(Am))' < OO. Combining this with we get Em>mi S'^l^-lAm) - 

0_(Am)) < CX). □ 

To reveal some further properties of the limit limm->.oo 5™Am, the following lemma is 
required, which is a generalization of formula (15. ip . 

Lemma 5.2. Let m > 2. a^n, Pm be two consecutive eigenvalues in with am < Pm- 
Then V/ G N, 

0^^(am) <0^^(/3m). (5.3) 

Proof. First we need to prove the following relation. 

PmM^lHam)) ~ (-l)'-Vm+i(0-(«m)), V/ G N. (5.4) 

In fact, when / > 3, using the Laplace theorem to expand the determinant qm+i{(f)^- {am)) 
according to the last (^ — 1) rows, we have 

qm+l{(j)- (am)) = g/(</>^^(«m))gm+l(</>-(ttm)) " / (</>- ^ ("m) )?«-! (</>- ("m) )'^ (</>- ("m) )gm ("m) • 

Since gm(ttm) = 0, we have 

qm+li<P- (Oim)) = g«(0- ("m))gm+l(0-(am))- 

This equality also holds for / = 2 by instead using an expansion along the last row of 
qm+2{4>-\am))- Hence for each / > 2, we always have gm+«(0- («m)) = g«(0- ("m))gm+i(0-( 
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Then from Lemma B in Appendix, we have qi{(t)'- {am)) > 0, hence qm+i{4'^- {(^m)) ~ 
Q'm+i(0-(ttm))- By the relation between pm+i{x) and qm+i{x), we can easily get f l5.4p . 

Now we prove fl5.3p . When 1 = 1, flS.Sp follows from (15. ip directly. In order to use the 
induction, assuming (15. 3p holds for /, we turn to prove 

Suppose am and (3m are the /c'th and {k + l)'th eigenvalues in Vm respectively. Recall 
that in Lemma 4.4, we have proved that Pm+i{(p~{am)) ~ (_x)™+'=-i. Combining this 
with (15. 4p . we have 

Pm+i+i{<P^-'''\am)) ~ (-1)-+'=+'-^ (5.5) 

On the other hand, if we denote am+i = (f)^-{am) and (3m+i = 4>^-{Pm), then it is easy 
to see that and (3m+i are the /c'th and {k + l)'th eigenvalues in respectively. 

Lemma 4.4 says that 

Pm.+l+li^4c^rn.+l)) ~ {-l)"^+^+'^~^ (5.6) 

and 

Pm+l+l{<P-{Pm+l)) ~ (-l)-+'+^ (5.7) 

Furthermore, if we denote Pm+i+i = (f)-'^^\Pm), then Pm+i+i is the only root of Pm+i+i{x) 
located between (f)_{am+i) and 0_(/3m+/), i.e., 

(j)-{am+l) < (3m+l+l < (p-i/^m+l)- (5.8) 

Noticing that from the induction assumption, we have (j)^l'^^\am) < (t>~{Pm+i) since 
Pm+i = 4>^-{Pm)- Moreover, (15. 5 p and (15. 7p say that there exists at least one root of 
Pm+i+i{x), denoted by /3m+i+i, between 0^"^^^(a„,) and (j)_{/3m+i), i-e., 

^^!:''\am) < (3*m+i+i < <PA(3^+i). (5.9) 

Since (j)-{am.+i) = (j)-{(j)^l\am)) < (p-'^^\am), we have 

0_ («„+/) < (j)^!^^\am) < ^m+l+l < 0-(/3m+/)- 

Combing this with (15. Sp . from the uniqueness of /3m+i+i, we have f3m+i+i = f^m+i+i- Hence 
substituting it into (15. 9p . we finally get (f)^l'^^\am) < Pm+i+i, i-e., (j)^l'^^\am) < (p-~^^\(3m), 
which is the desired result. □ 

The following is an application of Lemma 5.2. 

Lemma 5.3. Let mi > 2, ami,Pmi be two consecutive eigenvalues in Vm^ with 
ami < Prrii- {am}m>mi is an infinite sequence related by ttm+i = 0-(ttm.),Vm > mi; 
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{/3m}m>mi is an infinite sequence related by Pm+i = (p-{/3m),^rn > mi. Then Vm > mi, 
«m < Pm- Moreover, 

lim S'^am < lim S'"/?^,- 



Remark. In SQ case, this is a direct result since (p-ix) is a definite strictly increasing 
continuous function. 

Proof of Lemma 5.3. Let m > mi. Since am = 
we have 

< < = (3m (5.10) 

by Lemma 5.2. Hence Vm > mi, < 

Now we prove lim^^oo 5"'am < linim^oo 5"'/3m- 
Let m > mi. Then from f l5.10p . we have 

Hence /S^ — Om > 0^~"^^~"^''(0-(ami)) — 0^~"^^~"^''(0-(«mi))- Since > | whenever 

< X < 5, and < </)_(am.j) < 0-(ami) < 5, we have 

Hence 5'"(/9,„ — a^n) > 5™^+^(0_(amJ — which yields that 

hm 5"'(/3^ - a^) > 5'"^+i(0_(«„J - 0„(a™J) > 0. 

Thus hm„^oo 5'"am < hm^^oo 5'^/3m. □ 
Lemma 5.4. hmm^oo5'^Am > 0. 

Remark. In SQ case, this is a direct result, since {5'^Xm}m>mi is then a monotone 
increasing sequence. 

Proof of Lemma 5.4- Without loss of generality, we assume that Xmi is the least 
eigenvalue in V^^, since Lemma 5.3 says that it suffices to prove for this special case. 
Then Vm > mi, is also the least eigenvalue in V^. Note that Lemma B in Appendix 
says that Vm > mi, we have Am > 4>^T\q). Hence 

lim 5'"Am > lim 5'"0L"^(6) > 0, 

where the existence and positivity of the second limit are already shown in 5^ case. See 

m- ° 

Now we define 

A = ^ lim 5™Am. 

/ m— >oo 
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We will prove A is an primitive eigenvalue of A on the fractal fi. 

Note that Vm > uiq, G P+, i.e., is a root of both Pm{x) and qm{x) by Lemma 
4.5 and Theorem 4.1. As in Section 4, denote by A,-™'' the successor of Am, of order (m — i) 
with 2 < i < m for simplicity. Lemma 4.6 and Theorem 4.1 say that the system (14. 2 p 
of equations has 1-dimensional solutions {bi,b2,--- ,&m-i) with &i 7^ and fcm-i 7^ 0. 
We normalize the solution by requiring bi = 1, and write it as {b^i^\ b^\ ■ ■ ■ , b^^~i) with 
^(m) _ ]^ specify its relation to A^. We always denote feg™'' = for convenience. As 
described in Section 4, from {b^^\b^^\ ■ ■ ■ one can recover the unique (up to a 

constant) Am-eigenfunction Um on fi^ (noticing that X^"^^ 7^ 2 or 5, V2 < i < m). Hence 

UmlaQr,, = 0- 

For each m > rriQ, we start with the Am-eigenfunction Um on Q^, and extend Um. to Q 
by successively using the eigenfunction extension algorithm corresponding to the revised 
eigenvalue sequence {Am, 0-(Am,), (Am), ■ ■ ■ } (starting from Am, but continued with 
the standard spectral decimation eigenvalues) to get a primitive eigenfunction (possessing 
the symmetry in each cell Fyj{SQ) under the reflection symmetry fixing -F^^j^o with word 
w taking symbols only from {1, 2}) on Q. We still denote Um for this function. Obviously, 
Um may not satisfy the Dirichlet boundary condition. \/i > m, we use Xf^^ = 0^ ™''(Am) 
to denote the i-level revised eigenvalue. Hence for each m > mo, Um is an eigenfunction 
associated to the eigenvalue sequence {X[^^}i>2, where A^'"'' = (Am), V2 < i < m, 

and Aj-™'* = (j)^l~"^\Xm) , > m. We use bf^^ (Vi > m) to denote the value of Um at 
vertex F^qQ. Hence {bl"^^}i>o are the values of Um on the skeleton of Q which conversely 

K>o- 



determine Um on Q. We have the following relationship between {A^'"''}i>2 and {bf^^}i 



'U-SX^t^ + {6-)i: l)ti 
•'-^^^ (2-ASS)(5-AS) 



which follows from the eigenvalue equation at the vertex F^qo. Note that when 1 < i < 
m — 1, these are exactly the equations in (14. ip . Moreover, Um on Q satisfies that 



-AWm = 5"'<l>(Am)Wm OU Q, 

Um{qo) = 0, 

Um\L = hmj^oo^-™^ < 00, 



where ^{z) is a function defined by $(z) = | limfc^oo 5'^0^^(z). The existence of the limit 
limj^oo bf^^ will be given later. 

It is easy to find that 5™$ (Am) -H- A as m goes to infinity. For each m > mi, let 
Vm = II ""1*1 . We will prove that |fm.|m>mi contains a subsequence which converges 
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uniformly to a continuous function v on Q and w is a Dirichlet eigenfunction associated 

to A. We need the following lemmas. 

Lemma 5.5. There exists a constant Ci > depending only on mi, such that^i G N, 

G N, we have {b^^^ — < Ci(^)*||Mm||oo uniformly on m > mi. 

Proof. Without loss of generality, assume i > mi and A^i is not the largest eigenvalue 

in P^j. Denote by the next eigenvalue of A^i in V^^. Let {7m}m>mi be the infinite 

sequence staring from related by 7^+1 = 0-(7m), Vm > mi. We now show 

aS:;^? < 7^+1 < 0-(2), Vm > mi. (5.12) 
In fact if m > z + 1, then 

aH = f^-'-'-'KXm) = /('"-^-^H0L™"^"'Ha.+O) < /("-^-^H0^'""^"'^(A.+i)) = A.-,! < 7.+1. 

If m < i + 1, then A^-^] = (p^l^^ ™''(Am) < 0^^^ "^\lm) = li+i by using Lemma 5.2. The 
right inequality of fl5.12p is obvious. Hence ( I5.12p always holds. 
On the other hand, notice that from f l5.1ip . 

^(m) _ ^(m) ^ - 1,0 - \+i)Oi^i _ ^(m) 

' " 2(2-AS2)(5-AS:J) 



(6 - xt:l)it' - ti) - (20Aa - 9(Ai:;)^ + {x^ir)b^^ 
2(2-a1:;)(5-a;S) 



Hence 



(m) _ ,(m)| . \Q-Xi+l\ I (n^) _ {„) , |20 - 9X\'^l + (A,^!^'])^! (^n,) (^) 



In the remaining proof, we use c to denote different constants. 
By fl5rT2|) . we have 

(2-7i+i)(5-7i+i) 
Noticing that 7j = 0(i) and < ||Mm||oo, we get 

\i(Tn) i(m)i ^ / 3 '^Mjfm) i , II II 



Hence 



|,(m) _ ,(m)| ^ _ + il 



\U 



m oo- 



Similarly we have the estimates 



\i(m) i(m) I ^ <J li,(™) I I ^ II 
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till 

|,(m) _Am) I < ,{m)\ . |L, 11 



A routine argument shows that 



|,(m) ,(m)| , / 3 ,j_^j|,(m) , (m) i , / ^ 



j^gy I mi+i "mi I ' ^XO 5mi+l 

Hence we have proved that 



\Ur 



1 7 (m) 7 (m) I ^ / 3 N j II 



10' 



where c depends only on mi. 
Similarly, we have 



\bt:i-bt:i\<c{^r'\\uA^ 



till 



|,{m) _ ,{m) I ^ |'jiy+p-l|| II 



By adding up the above estimates, we finally get — l>(^\ < Ci(:^)'||M„J|oo-n 

Lemma 5.6. For each m > nii, limj^oo exists. Moreover, there exists a constant 
C2 > depending only on nii, such that | limj_^oo ^1'"'* I < C2(^)™||wm||oo uniformly on 
m > nil. 

Proof. For each m > mi. Lemma 5.5 says that each sequence {b["^^}i>i is a Cauchy 
sequence, hence limj_j.oo bf"^ exists. 

Taking z = m, p = 1 in Lemma 5.5, noticing that bm^^ = 0, we get that |fc^i| < 

r',l'A'|m|L, II 
^1 V 10 / II "mil 00 • 

On the other hand, Vi > m + 1, notice that \bf^'^\ < {b^/^'^ — b^m+il + l^m+il- By using 
Lemma 5.5 again, we have 

\bt^\ < Ci(^r+in^||oo + Cii^rWumWoo = C^2(^ri|n^||oo. 
Letting i — 00, we get the desired result. □ 

In the following context, let 6m denote the limit limj^oo ^I'^Vll'^mlloo- Lemma 5.6 

3 ^m J „, Um 



guarantees the existence of this limit, and furthermore, \6m\ < C*2(^)™'- Let w 
Vm > mi. Then Vm on Vt satisfies that 

-Ai;„ = 5"'$(Am)wm on Vt, 
Vm{qo) = 0, 

'^m I L ^m • 

Lemma 5.7. {dnVmi%)}m>mi is uniformly bounded, i.e., there exist a constant C3 > 
depending only on mi, such that \dnVm{(lo) \ < C3. 



43 



Proof. Let m > mi. Choosing a harmonic function h such that /i(go) = 1, h{Fiqo) = 
h{F2qo) = 0, the local Gauss-Green formula on Fq{SQ) says that 

^Fo{sg){vm,h) = - {Avm)hdiJ,+ ^ hdnVm. 

JPoisg) dFoisg) 

Hence \dnVm{qo)\ < \£Fo(Sg){Vra, h)\ + I Jp,^^^^g^{AVr,^)hdfl\. 

Since /i is harmonic on Fq{SQ), we have i^Fo{5g)(^^m5 h) = ^£{vmoFo, hoFo) = ^£o{vm.o 
Fo, hoFo). Noticing that h{qo) = 1, h{Fiqo) = h{F2qo) = 0, we get \£Fo{sg){vm, h)\ < Ci, 
since ||t;m||oo = 1- 

On the other hand, since Av^ = —5"^^{\m)vm, we have | Jp,^^gg^{Avm)hdfi\ < 5'"^^{Xm)\\'Vm\\oo- 
\\h\\^fi{Fo{Sg)) < C2, since 5™<I>(A„) ^ A. 
Hence |9„Wm(go)| < ci + C2 = C3. □ 

Lemma 5.8. {S{vm)}m>mi is uniformly hounded, i.e., there exists a constant > 
depending only on mi, such that S{vm) < C4. 

Proof. Wn > mi, let Kn be the part of Q above dQn\{<lo}- We first prove {SK„{vm)}m>mi 
is uniformly bounded and the upper bound is independent on n. 

Fix n > nil, ^ > f^i- The Gauss-Green formula says that jj^ Avmd^ = YldK ^nVm- 
From the symmetry property of f^, dnVm takes same value along dK^ \ {qo}- Hence we 
get 

-5'"$(A™)/ Vmd^l = ^nV„Mo)+2^^nVm{F^{qo)). (5.13) 



Kr 



On the other hand, the Gauss-Green formula also says that 

^K„{^m) = - {AVm)VmdfI + ^VmdnVm 

JKr,. 



= 5™$(A„) / vldfi + 2"vUFl'qo)dnVm{Fiqo), 
since Vm{qo) = 0. Combined with fl5.13p . it follows 

^i^„(Wm) =5'"$(A^) / vl^dfI + Vm{Fl'qo){-5"'^{Xm) V^d^l - dnVm{qo)) ■ 

Since 5™$(Am) — A, there exists a constant c > 0, such that 5''"$(Am) < c. Hence 

^KniVm) < c\\Vm\\lo + || f m || oo (c|| f m || oo + | C^n^ m (^o) | ) • 

Using Lemma 5.7, we get Sx^Vm) < 2c -|- C3 = C4. Since the above inequality is inde- 
pendent on n, we then get the desired result by passing n to infinity. □ 
Now we come to the main result of this section. 
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Theorem 5.1. There is a subsequence of {fm}m>mi; converging uniformly to a 
continuous function v on Q. Furthermore, v is a Dirichlet eigenf unction associated to 
A = I liiiim_j.oo 5'"Am- 

Proof. For each m > mi, since Vm G domS, we have 

\vUx) -v^{y)\ < \/x,y G Q, 

where d{-, ■) is the effective resistance metric on Q. Hence by Lemma 5.8, 

\vm{x)-vUy)\ < Cl'^d{x,yYl\ \/x,y G Q 

holds uniformly on m > mi. Thus {vm}m>mi is equicontinuous. Moreover, notice that 
{vm}m>mi IS also Uniformly bounded. Then using the Arzela-Ascoli theorem, there exists 
a subsequence {vmk} of {vm} which converges uniformly to a continuous function v on Q. 

Let Gfi{x,y) denote the Green's function associated to Q. See the explicit expression 
for Gn{x, y) in [12]. Then Vfc, we have 

/ Gn(x,i/)5™^$(A {y)dfi{y) + hn,,{x), (5.14) 

Jq 

where /i^fe is a harmonic function on Q taking the same boundary values as Vm^, i.e., 
hmkilo) = and h^^L = 6'^^.. If A; — )■ oo, then 9^,, — )■ and hence km,, goes to 
uniformly on Q by the maximum principle. Hence by letting /c — i- oo on both side of 
f l5.14p . we get 

v{x) = / Gn{x,y){Xv{y))dfi{y). 
Jn 

Thus we finally get 

—Av = Xv in f2, 

v\n = 0, 

i.e., i; is a Dirichlet eigenfunction associated to A. □ 

Thus for each sequence {Xm}m>mo^ ^6 have proved that A = | lim^^ 5"^Xm is a primi- 
tive Dirichlet symmetric eigenvalue of A. We denote by the totality of all this kind of 
eigenvalues. Lemma 5.2 and Lemma 5.3 guarantee that all these eigenvalues are distinct 
and they are all greater than 0. The skew-symmetric case is similar. We denote by V~ 
the set of skew-symmetric eigenvalues generated in this way. Let V = U V~ denote 
all the primitive Dirichlet eigenvalues of A. 

6 Complete Dirichlet spectrum of A 

It is clear that the weak spectral decimation recipe constructs many primitive eigenval- 
ues (hence also many miniaturized eigenvalues) of A. Recall that the standard spectral 
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decimation recipe also constructs many localized eigenvalues of A. It is natural to ask 
do these recipes construct all of the spectrum of A? In this section, we will answer this 
question. 

Till now, for each m > 2, we have proved that the spectrum Sm of the discrete 
Laplacian A^ consists of Cm, T^m and Aim the three types of eigenvalues. After passing 
the approximation to the limit, we have proved that there are at least three types of 
eigenvalues £, V and Ai in the spectrum S of A, i.e., S D C U V U JH. We call all of 
the above three types of eigenvalues raw eigenvalues. By the raw multiplicity of the raw 
eigenvalue A, we mean the multiplicity of the associated eigenvalue Xmo of A^o, where mo 
is the generation of birth. Since linear independent eigenf unctions of A^p belonging to 
Xmo Si^6 ^ise to linearly independent eigenfunctions of A, and the fact that all primitive 
graph eigenvalues have only raw multiplicity 1, the raw multiplicity of A is not greater 
than the true multiplicity of A. 

Denote by S' the collection of raw eigenvalues of A, then S' = CUVUAi and S' C S. 
Hence we need to prove S' = S and the raw multiplicity of each element of S' coincides 
with its true multiplicity. 

Comparing to the proof of the same problem in the SQ case, we only need to prove 
the following result. See details in [10]. Recall that = ]\{V^ \ dVtm) = ^^^^^ - 2"+^ 

Theorem 6.1. Let < ni < K2 < ■ ■ ■ be the rearrangement of elements of S' 
each repeated according to its raw multiplicity. Let {Km.,i}i<i<am be the m-level graph 
eigenvalues of Am. on Qm including multiplicities. Then 



In order to prove this theorem, we first list some notations and lemmas. It is more con- 
venient to consider the following slightly different classification of all the raw eigenvalues 



where V~ = V~ UAi, since miniaturized eigenvalues have the same generation mechanism 
as the skew-symmetric primitive eigenvalues. In the following, we always use a,/3,7 
to denote C, V~^, V~ type eigenvalues respectively. Accordingly, Vm > 2, all the m- 
level graph eigenvalues are classified into the three types Cm, T^m where "P" = 

Vm U J^m, and we always use am,(3m,lm to denote them respectively. For simphcity, 
we denote Am = %Cm, Bm = P^m and Cm = ^V'- Of course, am = Am + Bm + Cm- 
Moreover, recall that Pm{5) and p^(6) are the multiplicities of m-level initial eigenvalue 
5 and 6 respectively. See the exact values of them in Section 3. 




of A, 



s' = cuv^ uv 



Lemma 6.1. C = IJ^g (disjoint union) where C} 



C [5' 



.'=<l>(3),5'=$(5)]. 
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Proof. Va G £, let {am}m>mo be the corresponding sequence of eigenvalues with a 
generation of fixation nii. Then a = | linim^oo S'^Om = 5™i$(amJ. 

If a^ni is an initial eigenvalue, then can only be equal to 5. If continued 
eigenvalue, then ami — 0+('^mi-i), which yields that 3 < ami — 5- Hence we always have 
3 < ami < 5. 

Noticing that each localized eigenvalue has generation of birth at least 3, denote by 
the set of eigenvalues with mi = k, k = 3,4, Then £ = IJfcla C'' C 

[5^$(3),5*^<l>(5)]. Since 0_(5) < 3, we have $(5) < 5$(3). Hence £ = U^s'^'' is a 

disjoint union. □ 

Lemma 6.2. P+ = (disjoint union) where V^'^ C (0,5^$(6)] and V^''' C 

[5^$(0„(3)),5^<l>(6)] /or A; > 3. 

Proof. Wf3 G ■P"'', let {/3m.}m.>mo be the corresponding sequence of eigenvalues with a 
generation of fixation mi. Then /3 = | limm-^oo S^/Sm = 5™^ lim„^oo |5"</)^''(/3mi)- 

If (3mi is a continued eigenvalue (hence mi > 3), then we must have (3mi = 4>+{/3mi-i), 
which obviously yields that f3mi > where Pmi-i denotes the largest eigenvalue 

in Vmi^i- If Pmi is an initial eigenvalue with mi > 3, then obviously Pmi > 0-(/3mi-i)- 
Hence we always have /3mi > <P-{l^mi-i) if "^1 — 3- 

Moreover, When mi > 3, if we denote /S^i-i largest eigenvalue in Vmi^i except 
for then we have 0_(/3;^^„i) > 0_(/3;^*^_i). It is easy to check that > 

(p+{/3mi-2) > 3 since mi > 3. Thus f3mi > 0_(/3j^^_i) > 0-(3). When mi = 3, it can 
be checked directly that /^s > 0_(/32) ~ 1-33 > 0-(3). Hence we always have f3m.^ > 
4>~{l^mi-i) > 0-(3) if mi > 3. By Lemma 5.2, we have 0L"^(/3^J > 0L"^(0_(A;;,_i)), Vn G 
N. Hence if mi > 3, we have 

/3 = 5™i lim |5>L"^(/3„J 

n— >oo / 

> 5- lim |5>L")(0„(/3;,^„i)) 

> 5™^ lim -5"0L"^(0_(3)) 

n—>oo 2 

= 5'"^<l>(0_(3)). 
On the other hand, when mi > 2, we always have 

/3 = 5"'' lim -5"0L"H/3mJ < 5"^^ lim -5"0L"^(6) = 5'"i$(6). 

n— >oo 2 n— >oo 2 

Denote by the set of eigenvalues with mi = k, k = 2, 3, ■ ■ ■ . Then = 

U^2^"^''' where P+'^ C (0,52$(6)] and C [5'=<l'(0_(3)), 5'=<l>(6)] for k > 3. 

Next we need to prove = [J^=2'^^'^ is a disjoint union. \/2 < k < k', take an 
element /3 in 'P+''^, /?' in 'P+''^' respectively. Then /3 = S'"' lim„_j.oo |5"0^''(/3fc) for some 
eigenvalue Pk in "P^, and /?' = 5'^' lim„_^oo |5"</)^^(/3^/) for some eigenvalue in V^,. 
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Note that ^L*'' ^\(3k) and both belong to V^,. Since k' is the generation of fixation 
of /3', we can easily get 0^ '^\f3k) < P'k'- Then by using Lemma 5.3, we have (3 < f3'. 

From the arbitrariness of /3, /3' and k, k', we finally get that = UfcL2^^''^ ^ 
disjoint union. □ 

Lemma 6.3. Let < ai < 02 < ■ ■ ■ be the rearrangement of elements of C each 
repeated according to its raw multiplicity. Let {«m,i}i<i<yim be the m-level localized eigen- 
values of Am on Qrn. including multiplicities. Then 

lim 



1 

providing Y.7=\ i < ^• 



\<i<Ara 2 i=l * 



Proof. Noticing that lim,„_>.oo ^"gif^ = 0, it suffices to show that 

^m-p^(6) 



y — y (6.1) 



i<»<A™ 2"' i=l " 



converges to as m goes to infinity. 

Vm > 2, denote = — Pmi^)- By Lemma 6.1, {ai, 02, ■ ■ ■ , aD„i} is an arrange- 
ment of elements of IJfcLs each being repeated according to its raw multiplicity. The 
first sum of (16. ip has also Dm terms, which can be rearranged so that 

lim J5"^+"0L"^(«™,) = 1 < 2 < Dm. 

Hence by using Lemma 6.1, (161!) is equal to YlT=^^ry (^fkiisr:^ 77-)- If '^i ^ ^'^ 

(/c = 3, ■ ■ ■ ,m), then = 5'^<l>(^) for some 9 G [3,5] and accordingly the corresponding 
am,i is of the form am,i = (t)^^~^\0). Hence 

< 1 _ 1 . 1 1 . 

^ ~ ~ 5^ |5— '^(/.L™"*^^^) ~ *M 

Since 3g„^'(n)^- ^ converges to uniformly on [3, 5] as n goes to infinity, Ve > 0, the last 
expression is dominated by ^ whenever m — k is greater than some number A^. When 
m — k < N, the same expression is dominated by for R = | inf3<a.<5 The 
number of a^'s in is less than Ak^i + p^{5), so (16.11) is dominated by 



-e 



I^k^l^fj^ Om, 

^ 5^ - ^ ^ <cie + C2(-) - 

fe=3 k=m-N 

for some constants Ci,C2 > 0. Then let m be large enough, (16.11) can be dominated by 

(ci + 02)6. Hence we have proved y'i<i<A„ -^—^ Ylf-^ — converges to as m goes to 

infinity. □ 
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Lemma 6.4. Let < /3i < /32 < ■ ■ ■ be the elements of in increasing order. Let 
{Pm,i}i<i<B,n be the m-level primitive symmetric eigenvalues of on Q^- Then 



1 °° 1 

l<i<Bm 2 1=1 ^* 



providing Y^Zi jl < ^■ 

Proof. It suffices to prove that 

converges to as m goes to infinity. 

By Lemma 6.2, {/3i,/32, ■ ■ ■ ^/^b^} is an arrangement of elements of UfcL2'^^''^- "^^^ 
first sum of (16. 2p can be rearranged so tliat 

lim ^S'^+^L"^/?^,) = A, 1 < ^ < B^. 

Hence by using Lemma 6.2, (16. 2p is equal to Y.T=2 Z]ftg-p+,fc( 35^^ - j-)- The /c = 2 



term converges to as m goes to infinity since ^V'^''^ = B2 = 3. 
Hence we only need to prove 

converges to as m goes to infinity. 

If (3i E (A; = 3, ■ ■ ■ ,m), then (3^ = 5Mim„^oo |5"0^"V) for some 9 E V,^ and 

accordingly the corresponding Prn,i is of the form Prn,i = ''\6). Hence 

11,1,1 1 , , , 

(6.4) 



'§5-/5^,, A' S'^'lsm-fc^C— ^^)(^) lim„^oof5"0L"^(^) 
From the proof of Lemma 6.2, we have 

^5>L"H0-(3)) < ^5>L"V) < ^5>L")(6). 

Then by the proof of Lemma 5.1, Ve > 0, the right side of formula (16. 4p is dominated by 

4k£ whenever m — is greater than some number A^. When m — k < N , - — ... — is 

s*- — ' Is™;/.!:! '(6») 

dominated by ^;1^ for R = |0'l^^^''(3). The number of /3j's in 'P+''^ is controlled by fi^, 
so the sum (16. Sp is dominated by 

fc=3 k=m-N k=m~N fj^eV+'l' 

9 1"" 1 

s ^■- + '^(5)'"5+ EEs- 

k=m-N ftGP+''= 
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Noticing that X^i^i i" < ^^^^ term goes to as m goes to infinity. Hence for large 



/3. 

m, (16.51) is less than (ci + C2 + l)e. Thus we have proved X]i<i<B„ 35^^ - Z^i^ 
converges to as m goes to infinity. □ 

Lemma 6.5. Let < 71 < 72 < ■ ■ ■ be the elements ofV~ in increasing order repeated 
according to their raw multiplicities. Let {7m,j}i<i<Cm be the m-level V~ type eigenvalues 
of Am on Qrn including multiplicities. Then 



1 °° 1 

lim =y — , 



providing YlZi ^ < 

The proof is similar to those of Lemma 6.3 and Lemma 6.4. 

Proof of Theorem 6.1. Let < Ki < K2 < ■ ■ ■ be the rearrangement of elements of 
S each repeated according to its true multiplicity. Let Vi,V2, ■ ■ ■ be the associated eigen- 
functions. Let Gn{x,y) be the Green's function for Q. Then Gn{x,y) can be expanded 
uniformly convergence series 

1=1 

Since S' G S and the raw multiplicity is not greater than the true one, we get that 
< Hence < E»=i ^ < ESi ^ < The by adding up 

the results in Lemma 6.3, Lemma 6.4 and Lemma 6.5, we have 

1 °° 1 

lim y = > — < oo.n 

l<.l<am ^ «=1 

Based on Theorem 6.1, using similar argument in [10], we finally get S' = S and the 
raw multiplicity of each element of S' coincides with its true multiplicity. Thus we have 

Theorem 6.2. S = £ U P U where the union is disjoint. 

Hence we have constructed the complete Dirichlet spectrum of A on Q. 

Next we describe the Weyl's eigenvalue asymptotics on Q. We will find a formula 
analogous to (11. ip . We define the Dirichlet eigenvalue counting function 

p^(x) = l{\eS : X<x}, 

repeated according to multiplicities. Then we have: 

Theorem 6.3. There exist positive constant c, C such that 
for all X large enough, where ds = log9/log5 is the spectral dimension of SQ . Moreover, 

< liminf p^(a;)a;"°'^/^ < limsupp^(a;)a;"'^'^/^ < 00. 
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Proof. We divide p^{x) into four parts p^{x), p^^(x), (x) and p^{x) corresponding 
to different types of eigenvalues. The exact definitions are: p^{x) = jj{A G £ : A < x}, 
pP+[x) =^{\eV^ ■.\< x}, p^lx) = i{\eV~ : \<x} and p^{x) = ^{X E M : \ < 
x}. Obviously, 

= p^{x) + p^\x) + p^~ {x) + p^{x). 

For p^(x), it is same as the SQ case, hence there exist positive constant c', C such 
that dx'^^l'^ < p'~'{x) < C'x^'^^'^, for all x large enough, and furthermore 

< liminf p^{x)x~'^'^^'^ < limsup p^{x)x~'^^^'^ < oo. 

See details in [TU] . 

Next we consider p^^ (x). Denote the largest eigenvalue in and Z?^*"-* the eigen- 
value in P"*" corresponding to the sequence {4>^-\P^)}n>o, i-e., Z?'-™''' = lim„^oo |5"'^"^0^'*(/3^). 
By using Lemma 5.2, it is easy to check that 

015"^ = lim -5"+'^0L"H2) < lim -5"+'^0L"H < < lim -5"+'^0L"H6) = 025"^, 

n^oo 2 n— >-oo 2 n— 5>oo 2 

(6.6) 

for appropriate constants ci, C2 > 0, where denote the largest eigenvalue in except 

B* 

Notice that the bottom eigenvalues in are generated from eigenvalues in by 
extending these eigenvalues by choosing 0_ relation for all m' > m. Hence we get 

p^+(/3M)=r„, Vm>2. 

Using fl6.6p . we get p^^(ci5"^) < r^, and p^^(c25™) > r^. 

Denote by /cq the least number such that 5'^"ci > C2. Vx > 25c2, choose a number m 
such that C25™ <x< C25™+^ Then C25™ < x < ci5'"+'=o+^ Hence 

,3^1og2/log5 < < p^+(c25™) < p^\x) < p^"(Ci5-+'=«+^) < r^+,„+, < C4x'°^^/^°^^ 

for appropriate constants C3, C4 > 0. Thus we have proved that for x large enough. 
Similar argument yields that for x large enough, 

^^^log2/log5 < pP-(^^^ < ^^^log2/log5^ 

for appropriate constants c^.Cq > 0. 
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Now we consider p^{x). Notice that for each A' G {A G : A < x}, there exists a 
k >1, such that A' has multiphcity 2^ in A^, and ^A'g{Ag'P~:A<-^}. Hence 

p-(x)<$:2v-(j). 

Denote A* the least eigenvalue in V~ . Then 

[log(x/A, )/ log 5] [log(x/A. )/ log 5] 

k=l k=l 

for an appropriate constant Cy > 0. 

Taking the above estimates into account, we finally get the desired result. □ 



7 The Neumann case 

In this section, we give a brief discussion on the Neumann spectrum of A. It suffices 
to make clear all the discrete Neumann spectra of A^'s. As indicated in Section 3, we 
want to impose a Neumann condition on the graph Qrn by imagining that it is embedded 
in a larger graph by reflecting in each boundary vertex and imposing the Am-eigenvalue 
equation on the even extension of Um- It is convenient to allow m = 1, in which case 
there are only three boundary points in Qi and no others. Denote by the totality 
of primitive Neumann eigenvalues of the discrete Laplacian A.^. Due to the eigenspace 
dimensional counting argument in Section 3, this time we need to find out 2"* symmetric 
primitive eigenvalues and 2™ — 1 skew-symmetric primitive eigenvalues. 

We focus our discussion on V^'^ , the symmetric case, and describe a similar weak 
spectral decimation which relates "Pm'^ with Pm+i- Let Um be a Am'GigGnfunction of A^ 
with A„ G Still denote by {bo, bii ' ' ' , the values of Um on the skeleton of 

Write aS""^ the successor of A^ of order (m — i) with 1 < i < m. (This time we begin 
with X^^\) Assume that none of Xf^^^s is equal to 2 or 5 for 2 < ? < m. Then Um is 
uniquely determined by {bo,bi,--- ,bm). In addition to the eigenvalue equations at the 
vertex Fiq^, F^qq, • • ■ , F™~^go as described in Section 4, we impose the equations 

(4-AS'"))6o = 46i (7.1) 

at go and 

(4 - XJbm = 2bm-l + 2bm (7.2) 

at FI^Qq according to the Neumann boundary condition. Hence {bo,bi,--- ,bm) can be 
viewed as a non-zero vector solution of a system of equations consisting of m + 1 equations 
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in m + 1 unknowns, whose determinant is 



(m) 



(m)^ 



s{\t^) r(AL™)) 



-2 2 - A 



(m) 



Hence should be a solution of the following equation 

4_/(™-i)(a;) _4 

s(/("-2)(a;)) r(/(™-2)(x)) 



/(x) s(x) r(x) 

-2 2-x 



(7.3) 



Thus if Am is a root of q^^{x) and none of /^*^(Am)'s with 0<z<m — 2is equal to 
2 or 5, then A^ G V^^'^ . We should mention here that when m > 2, comparing to ^^(2^) 
in the Dirichlet case, q^{x) is a (m + 1) x (m + 1) tridiagonal determinant, containing 
qm{x) in the center as a (m — 1) x (m — 1) minor, i.e., we have 



4_/(™-i)(a;) 






-4 




















r(x) 
-2 2-x 



The degree of qZ{x) is 3(2™"! - 1) + 2'"-^ + 1 = 2™+i - 2, since the degree of qm{x) 
is 3(2™"^ — 1). The analysis on q^{x) is more complicated than that on qm{x) since for 
qm{x) we can always use the expansion of qm{x) along the first or last row to get a relation 
with two polynomials in same type but with smaller degree. 

The following lemma is a slight modification of the form of q^{x) from a (m + 1) x 
(m + 1) determinant to a m x m determinant. 

Lemma 7.1. Let m > 2. Then 

4_/Cm-i)(,x) _4 

r(/(™-2)(x)) 



g^(x) = (2-x)(x-6) 



/(/(x)) sifix)) r{f{x)) 
1 /(x) - 1 
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Proof. Substituting the expression for r{x) into fl7.3p . we get 



N 



(x) = (2 - x) 



4-/(™-^)(x) -4 
/(/("^-2)(x)) r(/('"-2)(x)) 



-2(5 - x) 
-2 1 



(2-x) 



(2-x) 



/(/("-2)(x)) r(/('"-2)(x)) 

l{x) s(x)-4(5-x) -2(5 -x) 
1 

4_ f{rn-l)(^x) -4 
/(j(m-2)^^)) s(/(™-2)(x)) r(/('"-2)(x)) 

l{x) s{x) - 4(5 - x) 



Noticing that s(x) — 4(5 — x) = (x — 6)(/(x) — 1) and /(x) = x — 6, we get the desired 
result. □ 

The following lemma focuses on the possibility of the roots of q^{x) satisfying /^'^(x) = 
2 or 5 for some < i < m — 2. 

Lemma 7.2. Letm > 2, andx he a predecessor of 2 or 5 of order i withO <i< m — 2. 
Then q'^(x) = 0. 

Proof. Firstly, let x be a predecessor of 2 of order i with < i < m — 2. Then 
/W(x) = 2 and /(^+^)(x) = 6. If < i < m — 2, the proof is the same as that of Lemma 
4.2. So we only need to check the i = m — 2 case. In this case we have 



/(/(™-2)(x)) s(/('"-2)(x)) r(/('"-2)(x)) 



-2 -4 
-4 -8 



Secondly, let x be a predecessor of 5 of order i with < i < m — 2. Then f^^\x) = 5 
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and 



x] 



j{m — Q_ Hence we have 



4 -4 
/(O) s(0) r(0) 



/(O) s(0) r(0) 

/(5) s(5) r(5) 



4 -4 
-6 26 -20 



-6 26 -20 

-1 1 



Thus we always have g, 



N I 



X] 



0. □ 



This lemma means that for m > 2, all the predecessors of 2 or 5 of order i with 



< z < m — 2 are unwanted roots of q^{x 



9m (a;) 



To exclude them out, we define 

for m > 2, 



(x - 2)(x - 5) ■ ■ ■ - 2)(/(™-2)(x) - 5) 



and 



Now we can say if is a root of the polynomial 



N, 



X) 



then Am G V^''^ . It is easy to 



check that the degree of Pmi^) is 2™, since the degree of q^i.^) is 2^^^ — 2 and the number 
of all the unwanted roots of is 2(1 + 2 + ■ ■ ■ 2 



m-2^ 



2 for m > 2 and for 



m 



1. The following is an easy observation on 



Lemma 7.3. For m > 1, p^{x) always has roots and 6. 

Proof. We only need to check g^(0) = ^^(6) = 0. It is easy to see that 



9^(0) 



4 -4 

/(O) s(0) r(0) 



/(O) s(0) r(0) 
-2 2 



4 -4 
-6 26 



-20 



26 
-2 



-20 
2 



9m(6) = follows from Lemma 7.1 for m > 2, and from direct computation for m = 1. 



□ 



In order to study the distribution of roots of p^{x), we now introduce a type of 
auxiliary polynomials lm{x) associated to First, Vm > 1, let lm{x) denote the 

m X m minor located in the upper left corner of q^{x), i.e., /i(a;) = 4 — a; and for m > 2, 



4-/("^-i)(x) -4 

r(/("*-2)(x)) 



55 



Note that the (m — 1) x (m — 1) minor located in the bottom right corner of Imix) is 
qm{x). The degree of lm{x) is 2™+^ — 3 since it is reduced by 1 comparing to the degree 
of q^{x). With similar argument in the proof of Lemma 7.2, all the predecessors of 2 or 
5 of order i with < i < m — 2 are roots of lm{x)- To exclude them out, we define 

"^^^^ ~ (x - 2)(x - 5) ■ ■ ■ (/('»-2)(x) - 2)(/{'»-2)(x) - 5) ' - ^' 

and 

li{x) = li{x). 

It is easy to check that the degree of lm{x) is 2™ — 1, since the degree of lm{x) is 2™+^ — 3 
and the number of all the unwanted roots of Imix) is 2(1 + 2 + ■ ■ ■ 2™~2^ _ 2™ _ 2 for 
m > 2 and for m = 1. 
Based on the property 

Ux) = s(x)L-i(/(x)) -r(/(x))/(x)L-2(/('n^)), 

lm{x) can be analyzed in a similar way that of Pm{x) or Pmix) in the Dirichlet case. We 
then have: 

Lemma 7.4. /,„(0) > and /m(6) < 0, Mm > 1. 

Proof. /m(0) > follows from a similar argument in the proof of Proposition 4.1 and 
Proposition 4.2. 

To prove /m(6) < 0, we only need to prove /m(6) < by the definition of lm{x). It 
can be checked that /i(6) = —2 < and ^2(6) = —40 < by an easy computation. For 
m > 3, an expansion of /m(6) along the first row yields that 

L(6) = (4 - f^"'~'\6))qU6) + 4(/(— 2)(6) - 6)g„_i(6). 

Recall that in the proof of Proposition 4.1(3), we have proved that qm{Q) < qm~i{Q) < 0, 
Vm > 3. Hence 

L(6) < (4-/("'"i)(6) + 4/(™"2)(6)-24)g„(6) = (f(^~^\6) -5){f^"'''\6) +4)qmi6) < 0, 

noticing that f^'^^'^^G) < —6 whenever m > 3. □ 

Lemma 7.5. For each m > 1, lm{x) has 2*^ — 1 distinct real roots between and 6 
satisfying 

< /3m,l < /3m,2 < ■ ■ ■ < /3m,2'"-l < 6. 
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Moreover, 

< (3ra+l,l < 0-(/3m,l), 

< Pm+i,k<(l>4Pm,k), V2<A;<2'"-1, 

0-(/3m,2'"-l) < /3m+l,2'" < 0+ (/3m,2'"-l ) 5 
0+(/3m,2'"+i-fc) < /3m+l,fe < 0+(/3m,2'"+i-/c-l), + 1 < < 2'""'""^ — 2, 

0+(/3m,l) < /3m+l,2'»+i-l < 6. 

Proof. It follows from a similar argument in the proof of Lemma 4.4. □ 
The following lemma shows a relation between p^(x)'s and lm{x)^s. 
Lemma 7.6. Let m>2. Then p^{x) = (2 - - 4/„„i(/(x)). 

Proof. This is easy to get since we have 

qZix) = (2 - x)L(x) - 4(2 - x)(5 - Vm > 2, 

using the expansion along the last row of q^i^)- ^ 
Now we consider the distribution of roots of p^{x). 

Lemma 7.7. For each m > 1, Pmi^) has 2*" distinct roots between and 6 (including 
and 6). Moreover, Pm(0+) < and p^(6— ) < 0. 
Proof. When m = 1, it naturally holds. 

Let m > 2. From Lemma 7.5, lm{x) has 2"* — 1 distinct real roots between and 6 
satisfying 

< i3m,l < (3m,2 < ■ ■ ■ < /3m,2'"-l < 6. 

For each 1 < A; < 2™ — 1, using Lemma 7.6, we have 

P^Wm,k) = (2 - l3m,k)lmWm,k) " 4/„_i (/ ) = -4/^-1 ) • 

When /c = 1, by Lemma 7.5, < f3m^i < (p~{f3m-i,i), hence < /(/3m,i) < /3m-i,i- 
Combined with /m-i(O) > from Lemma 7.4, it follows lm-i{f {(3m,i)) > 0, hence p^(/3m,i) < 
0. 

When 2 < k < — 1, following from Lemma 7.5, we have <^_(/3m-i,fc-i) < /3m,fc < 
<^_(/3m-i,fc), hence P>m-\,k~\ < f{l^rn,k) < /3m-i,fc- Combined with /m_i(0) > 0, it follows 

lrn-lU{Prn,k)) ~ {-if-', heUCe ~ (-1)^ 

When k = 2™~^, following from Lemma 7.5, we have 0_(/3m-i,2'"-i-i) < /9m,2'"-i < 
2^-1-1)5 hence /(/3m,2'"-i) > /3m-i,2'"-i-i- Combined with /m_i(0) > 0, it follows 
(/(/3m,2™-0) < 0, hence p^(/3m,2— > 0. 
When 2"^^^ + 1 < /c < 2™ — 2, following from Lemma 7.5, we have 0_|_(/3m-i,2'"-fe) < 
(3m,k < 0+(/3m-i,2™-/c-i), hence /3m~i,2^~k^i < f{f3m,k) < /3m-i,2'"-fe- Combined with 
/„_i(0) > 0, it follows ~ hence p^iPm,k) ~ (-!)'• 
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When k = 2"^ — 1, following from Lemma 7.5, we have 0+(/3rrt-i,i) < /3m,2'"-i < 6, 
hence /(/3m,2™-i) < Pm-i,i- Combined with /m_i(0) > 0, it follows lm-i{.f{Pm,2"^-i)) > 0, 
hence p^(/3m,2™-i) < 0. 

Hence we have proved p^{(3 m k ) 

(-1)'=,V1 < A; < 2"^ - 1. So there exist at least 
2m _ 2 ];oots of each located strictly between each two consecutive /3m,fc's. Moreover, 

Lemma 7.3 says that and 6 are also roots of p^{x). Thus we have found 2"^ distinct 
roots of Pmi^)- Since the order of p^i^) is also 2™, these are the totality of roots of 

Furthermore, from the fact that < and p^(/3m,2'"-i) < 0, we havep^(0+) < 

and p^(6-) < 0. □ 

In all that follows, we denote 

Am.,1 = < Am,2 < ■ ■ ■ < Am,2'"-l < Am,2™ = 6 

the 2™ distinct roots of p^i^) in increasing order, Vm > 1. In order to study the relation 
of roots of two consecutive p^(x)'s, we prove the following two lemmas: 
Lemma 7.8. Let m > 1 and 2 < k < 2™, then 

N (, (\ \\ ^m,k(p-{^m,k) — ^ , N 
Pm.+l\'P-\'^rn,k)) — ~ _i fx ^ F ' 'ml^»n,fcj 



and 



2 0-(Am,fc)-5 

Pm+l\.(P+\.^m,k)) — ^"TT^ \ ? ■ ''mK.^ 



2 (j)+{Xni,k) - 5 



m,k) 



Proof. For simplicity we only prove the first equality. The second will follow from a 
similar argument. It is easy to see that A^,*; is also a root of q^{x) and none of /^*^(Am,fc)'s 
(0 < i < m — 2) is equal to 2 or 5. 

By Lemma 7.1, q^^^{(l)_{\m,k)) = (2 - (l)-{\m,k)){(p-iKi,k) - 6) ■ A where 

4 - -4 

/(/('"-'HAm,.)) s(/("-')(A^,fc)) r(/(™-2)(A„,fc)) 



/(A, 



m,fc J 



■5(Am,fc) 
1 



r'(Am,,fc) 
Am,fc 1 



Noticing that from q^{Xm,k) = 0, we have 

4-/(™"^)(A„,fc) -4 
/(/('""'nAm.fc)) ^(/('"-'nAm.fc)) r(/(™-2)(A„,,)) 

^(Am,fc) 



s{Xm,k) r{Xm,k) 
— 1 1 — \m,k/'2 



0. 
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The summation of the above two determinants yields that A = -^Y^lm{^m,k)- Hence 

which yields the desired result. □ 

Lemma 7.9. Let m > 1. Then (-l)^"^/m(A„,fc) > 0, VI < A; < 2"^. 

Proof. Let /3m, i, f3m,2, ■ ■ ■ /3m,2™-i denote the 2"^ — 1 distinct roots of lm{x) in increasing 
order as described in Lemma 7.5. Then by the proof of Lemma 7.7, we have 

Am,l = < /3m,l < Am,2 < /3m,2 < " " " < Am,2'"-1 < /3m,2™-l < Am,2™ = 6. 

Combined with the fact /m(Am,i) = ^m(O) > by Lemma 7.4, it follows the desired result. 
□ 

Now we can prove the following Neumann analog of Lemma 4.5, 
Lemma 7.10. For each m > 1, V^'^ consists of at least 2"^ distinct eigenvalues 
satisfying 

Am,l = < Am, 2 < ■ ■ ■ < Am,2™-1 < Am, 2'" = 6. 

Moreover, 

0_(Am,/t-l) < Am+l,fc < 0-(Am,fc), V2 < < 2^", 
0+(Am,2™+i-fc+l) < Am+l,fc < 0+(Am,2'"+i-fc)) V2'" + I < k < T^^^ — 2, 

0+(Am,2) < Am+1,2^+1-1 < 6. (7.4) 



Proof. Noticing that each root of p^{x) belongs to V^'^ , we only need to prove the 
results for the roots oip^{x). The first statement follows from Lemma 7.7. We now prove 
the second statement. From Lemma 7.8 and Lemma 7.9, we have 

P^^ii<PAK..k)) = -V t'n"''!"^ ■ ^-(^-.^) ~ -U>^rn,k) - (-l)^ V2 < A; < 2™, 

2 (/)_(Am,fc)-6 

and similarly, 

P^+l(0+(A,n,.)) = -%^t4M^ ■ ^-(^--^) ~ -UKn,k) - (-l)^ V2 < < 2™. 

Following the above facts and Lemma 7.7, we can list the signs of the values of Pm+i(^) 
at different point x in the following table. 



X: 0+ </)_(Am,2) </>-(Am,3) ■■• </'-(A,Ti,2™) 0+(Am,2"O •■■ 4'+{^m,3) 4'+i^m,2) 6- 6 

P^+lW: - + - •■• + + ••• - + - 
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Hence there exist at least 2™'^^ distinct roots of satisfying (17.41) . Moreover, these 

are the totahty of the roots of p^_^_i{x) since the degree of p^_^_l{x) is also 2™+^. Hence 
we get the desired distribution of roots of p^^_^_i{x). □ 

The estimate 4>+{Xm,2) < ^m+i,2"^+^~i < 6 in Lemma 7.10 can be refined into 

0+(Am,2) < Am+l,2'"+l-l < 5 (7-5) 

by using the following lemma. 

Lemma 7.11. Form > 2, let Xm,i = 0, Am,,2, ■ " " , Am,2™-i, Am,2'" = 6 be the 2'" 
distinct roots of p^{x) in increasing order. Then 

^m,k + Am^2'"-fc+l = 5, V2 < < 2"* — 1. 

Proof. From Lemma 7.1, it is easy to see that if q^^{x) = and x 7^ 2 or 6, then 
g^(5 — x) = 0. Obviously, each X^^k (2 < /c < 2"* — 1) satisfies this property. □ 
Now we come to the main result of this section: 

Theorem 7.1. For each m > 1, V^'^ consists of at least 2™ distinct eigenvalues 
satisfying 

Am,l = < \m,2 < ■ ■ ■ < Am,2'"-1 < 5 < Am,2™ = 6. 

Moreover, 

0-(Am,,A;-l) < Xm+l,k < 0-(A™,fc), V2 < A; < 2"", 

0+(A™,2™+i_A.+i) < A„+i,fc < 0+(A„,2'"+i-fc), V2'" + 1 < A; < 2'"+i - 1. (7.6) 

Proof. It follows from Lemma 7.10 and Lemma 7.11. □ 
The following is a Neumann analog of Lemma 4.6. 

Lemma 7.12. Let \„i be a root of p^{x), Um a primitive Xm-eigenfunction on VLm, 
and [bo, fci, ■ ■ ■ , bm) the values of u„i on the skeleton of Q^- Then b^ ^ and b^ 7^ 0. 

Proof. Without loss of generality, assume m > 3. We still use A-™'' to denote the 
successor of A^ of order (m — i) with 1 < i < m. From the definition of p^{x), A-™'' 7^ 2 or 
5, for each 2 < i < m. Vector (60, bi, - ■ ■ , bm) can be viewed as a non-zero vector solution of 
system (14. 2 p of equations and in addition the two Neumann boundary eigenvalue equations 
dm) and (O. 

Suppose bm = 0. Then from (17.21) . bm-i = 0. It is easy to check that the deter- 
minant of the remaining equations in m — 1 unknowns (60, ■ ■ ■ ,&m-2) is /m-i(A^i). 
Since (60, - ,fcm-2) should be a non-zero vector, we have lm~ii^lnli) = 0, hence 
^m~i{X^m-i) = 0- Noticing that from Lemma 7.6, we have Pmi^m) = (2 — Xm)lm.{Xm) — 
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4/m,_i(A^i). Hence we get that lm{^m) = since both Zm-i(A^i) and Pmi^m) are equal 
to 0. But this is impossible, since Lemma 7.5 says that if /m-i(A^i) = then /m(Am) 
could not equal to 0. Hence bm 0. 

On the other hand, if 6o = 0, then by substituting it into the system, noticing that 
none of Xf^^'s is equal to 2 or 5, we can get bi = 0, ■ ■ ■ , 6^ = successively, which 
contradicts to 6^ 7^ 0. Hence b^ ^ 0. □ 

This is the whole story of the symmetric case. The skew-symmetric case is slightly 
different but very similar. The result is shown in Section 3, but the proof is omitted. 

With similar argument in the Dirichlet case, we finally get 

Theorem 7.2. For each m >1, consists of 2"^ distinct eigenvalues satisfying 

^m,l = < Xm,2 < ■ ■ ■ < Am,2™-1 < 5 < Xm,2"'- = 6. 

A relation between V^'^ and V^^^ is shown in (17. 6p . Similar properties hold for V~'^ 
with 2"* replaced by — 1, and Xra,i > in that case. Each Xm G has multiplicity 1. 
Moreover, the Neumann spectrum of Am on Qm satisfies 

where the union is disjoint. 

We should mention that Lemma 7.12 and its skew-symmetric analog show that there is 
no primitive eigenfunction (or miniaturized eigenfunction) that is simultaneously Dirichlet 
and Neumann {D — N). Hence the only possible D — N eigenfunctions are localized 
eigenfunctions. This is same as the SQ case. 



8 Further discussion 

In this section, we discuss to what extent our method can be extended to other domains 
in SQ. In particular, we will focus on fi^. (0 < x < 1). It seems that we can analyze the 
spectrum of A on case by case following the similar recipe for the Qi case. However, it 
is hardly to develop a general method which is suitable for all cases, although we believe 
that we are clear about the structures of the spectra. We let denote the bottom 
boundary of Q^- Thus will be a Cantor set for generic x, and a union of intervals if x 
is a dyadic rational. We may assume without loss of generality that 1 < x < 1, for if not 
we may first solve the problem for ^22;, and then simply dilate the solution to Qx- 

For simplicity, we only discuss the Dirichlet spectrum of A. Obviously, it will suffice 
to describe the discrete Dirichlet spectra of A^'s for all m. Hence the first problem is 
how to define the graph approximations. Similar to Qi, the fractal domain Qx can be 
realized as the limit of a sequence of graphs Qx.m- More precisely, Vm > 1, let V^'' 
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be a subset of Vm with all vertices lying along or under removed. Let ^lx,m be the 
subgraph of restricted to V^'^. Denote by dflx,m the boundary of the finite graph 
Qx,m- It is easy to find that V^'^ \ dQx,m, dQx,m approximate to Qx and dQx as m goes 
to infinity respectively. See Fig. 8.1 and Fig. 8.2 for Qx and Qx,m. where x = 3/4. On 
flx,m the Dirichlet Am-eigenvalue equations consists of exactly tKKn'" \ dflx,m) equations 
in tl(Kri'^ \ 9^x,m) unknowns. We denote by Sm{x) the spectrum of A.^ on flx,m for each 
m > 1. Accordingly, Sm{x) should consists of (at least) three types of eigenvalues, denoted 
by Cm{x), Vrn{x) and A^m(x) respectively. Vm{x) can also be split into symmetric part 
V^{x) and skew-symmetric part We omit the precise definitions since they are 

obvious. To ensure that there is no other eigenvalue in iSm(x), the following eigenspace 
dimensional counting formula is hoped to be held. 




Fig. 8.1. 03/4- 



We now focus on a particular example ^3/4 to illustrate how to extend the recipe 
for VLi. We should be particular interested in the primitive eigenvalues. We begin with 
P,^(3/4), the symmetric case. It is convenient to define the skeleton of f2m(3/4) by 
(go, FiQq, FioFigo, ■ ■ ■ , FioF^-'^qo) for m > 3 and (go, i^igo) for m = 1 or 2. Let be a 
Am-eigenfunction of with Am G P^(3/4). Denote by (60, ^2, ■ ■ ■ , bm) the values of Um 
on the skeleton of ^3/4,^ where 60 = = by the Dirichlet boundary condition. It is easy 
to observe that when i > 2, the eigenvalue equation at the vertex Fio-Fi~^go is exactly same 
as that of Qi case with suitably reindexed. Hence the generation mechanism of primitive 
symmetric eigenvalues is quite similar to the Qi case. Based on this observation, one can 
easily find that tj'P+(3/4) = 2"^ — 2 for m > 2 by still using the weak spectral decimation 
method. A similar argument yields that jjP~(3/4) = 2™ — 2"^~^ — 2 for m > 2. 
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^3/4,1 



■3/4,2 



^3/4,3 



Fig. 8.2. The first 3 graphs, ^^3/4,1, ^^3/4,21 ^3/4.3 i^i the approximation to ^^3/4 with inside 
points and boundary points represented by dots and circles respectively. 



To verify the eigenspace dimensional counting formula, we only look at the first 4 
levels of approximations since the continued process is similar. 

When m = 1, the result is trivial since there is no inside point in 1^3/4,1. Hence 

tl5i (3/4) =0 = ri'''''*\ 5^^3/4,1. 

When m = 2, it is easy to check that there are only primitive eigenvalues. Hence 
tl52(3/4) = ^2^(3/4) + ^2^(3/4) = 2 + 1 = \ d^W- 

When m = 3, it is easy to check that there are 4 initial localized eigenvalues, i.e., 5 
with multiplicity 1 and 6 with multiplicity 3; there are 6 primitive symmetric eigenvalues 
and 4 primitive skew-symmetric eigenvalues; there is no miniaturized eigenvalues. Hence 
^53(3/4) = tt/:3(3/4) + P3'"(3/4) + ^3^(3/4) = 4 + 6 + 4 = ^V^'^' \ 9^3/4,3- 

When m = 4, it is easy to check that besides 1 ■ 2 + 3 ■ 1 = 5 continued localized 
eigenvalues, there are 18 initial localized eigenvalues, i.e., 5 with multiplicity 4 and 6 
with multiplicity 14. Hence |:j£4(3/4) = 5 + 18 = 23. There are 14 primitive symmetric 
eigenvalues and 10 primitive skew-symmetric eigenvalues. Hence jlP4(3/4) = 14 + 10 = 24. 
Moreover, there are some miniaturized eigenvalues which come from the miniaturizations 
of eigenvalues in V:^{1). Hence tl-M4(3/4) = 2 ■ V^il) =2-2 = 4. Thus tl'54(3/4) = 

It is easy to verify the general formula for general m. We will not attempt to list the 
details here. However, a more important fact should be pointed out is that for ^3/4 case, 
the miniaturized eigenvalues in A^m(3/4) are generated not from those in 7^^(3/4) but 
from those in V^{1) for k < m — 2. This means to study iSm(3/4), one should first make 
clear i5m(l). Things will be more complicated for general Qx- 

Next we briefly present another observation. Still consider a domain with a series 
of graph approximations {flx,m}- Notice that there are only two possible patterns when 
passing from the m-level graph approximation to its next level. One is that the boundary 
dflx,m+i remains unchanged, i.e., dQx,m+i = dQx,m, the other is that dQx,m\{<lo} becomes 
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a collection of inside points of Qx,m+i, i-e., each point in dQx,m \ {qo} is connected with 
two new (m + l)-level points in dQx,m+i- In fact, for the SQ case, when passing from one 
level to the next level, the boundaries of graphs are always Vq, keeping unchanged. This 
is also the reason why spectral decimation can work for 2-series eigenvalues (which should 
be considered as the primitive eigenvalues in SQ case). As for the Qi case, when passing 
from one level to the next level, the boundaries always change. Due to this phenomenon, 
the spectral decimation recipe should be replaced by the weak spectral decimation recipe 
for primitive or miniaturized eigenvalues since their supports always touch the boundaries. 
For general (0 < x < 1), these two possible patterns can both exist. It is natural to 
expect that under the first pattern, the two levels of primitive eigenvalues are related by 
the spectral decimation (it is obviously true.), while under the second pattern, they are 
related by a weak spectral decimation instead. Thus we expect: 

Conjecture 8.1. For a domain (0 < x < 1) with a series of graph approximations 
{^x,m}, if the boundaries change when passing from m-level to (m + l)-level, then there 
is a weak spectral decimation relating the two levels of primitive symmetric (or skew- 
symmetric) eigenvalues. 



9 Appendix 

Theorem A. For each m > 2, let Xm,i, Am,2, ■ ■ ■ , Xm,rm ^he r^, distinct eigenvalues 
in in increasing order. Then Am+i,r™+i > 2. 

To prove this theorem, we need the following lemma: 
Lemma A. ^2(2) < 0, ^3(2) > and (-1)>„(2) > 0, Vm > 4. 
Proof. It is easy to check that ^2(2) = —8 < and ^3(2) = 68 > 0. 
Let m > 4. Then 

Pm{x) -- 



(x-2)(/(a:)-2)---(/(— 3)(x)-2) 
s(/(""^)(a:))g^„i(x) - lif^"'-^\x))r{f^^~^\x))q^^2{x) 

(x-2)(/(x)-2)---(/(™-3)(x)-2) 
s(/(™-^)(x)) -/(/("'-3)(3,))^(;(m-2)(^)) 

-Pm-l{x) + ,,,,,,, oX/W^^:.W^N 7;^Pm-2{x). 



/(m-3)(a;) _2^"'-^^ ^ (/(™-4)(^) _2)(/(™"3)(^) _2)^ 

Noticing that /(/(™-3)(x)) = /('^-3)(x) - 6 = (/(™-^)(x) - 2) (3 - /(™-4)(x)) and 
choosing x = 2, we have 

.(/(-"^)(2))p^„i(2) + 2(2 - /(-^)(2))(5 - /(-^)(2))(3 - /(-^)(2))p^„2(2) 

(9.1) 
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We will prove the following stronger result than that stated in Lemma A. 

p^(2) ~ (-1)™ and p^+i(2) +p„,(2) ~ (-l)'"+\Vm > 4. (9.2) 

Using ([9IID, it is easy to check that ^4(2) = 14064 > and ^5(2) = -593514756 < 
by a direct computation. Hence (19. 2 p holds for m = 4. In order to use the induction, we 
assume ( 19. 2p holds for m and will prove it for m + 1. 

First, it is easy to get that Pm+i(2) ~ (—1)™+"^, since otherwise pm+i(2) + Pm(2) ~ 
(— 1)™, which contradicts to the induction assumption. Hence we only need to prove 

P^m+2i2) + Pm+li2) ~ (-1)-. 

Note that from ( 191|) . 

Pm+2{2) +Pm.+l(2) 

(5(/(-)(2)) + /(-^(2) - 2K+i(2) + 2(2 - /("')(2))(5 - /(-)(2))(3 - /(-^)(2)K(2) 

/{»n-l)(2) - 2 

(2) + 6m(|}„+i(2) + Pm(2)), 

where 

_ s(/("^)(2)) + /(™"i)(2) - 2 - 2(2 - /('^)(2))(5 - /(™)(2))(3 - f "'~^^{2)) 

- /(m-l)(2) -2 

and 

_ 2(2 - /('")(2))(5 - /(™)(2))(3 - /("^-^)(2)) 
~ /(m-i)(2) -2 ■ 

It is easy to check that 6„ < 0, since f^"'\2) < /('"-i)(2) < /('"-2)(2) < noticing 
that /(2)(2) = -6 and m > 4. We will prove that < also. In fact, the numerator of 
am is s(7) + /(/3) - 2 - 2(2 - 7)(5 - 7)(3 - where 7 = /("^)(2) and (3 = /('"-2)(2) for 
simplicity. By using 7 < /(/3) < /? < —6, it is easy to get 

s(7) + /(/3)-2 = (2-7)(4-7)(5-7)-14 + 37 + /(/3)-2 

> (2-7)(4-7)(5-7)-16 + 47 

> (2-7)(4-7)(5-7)-(2-7)(5-7) 
= (3-7)(2-7)(5-7), 

and 

3 - 7 > 3 - f{/3) = 3 - /3(5 - /3) > 3 - 5/3 > 2(3 - 

Hence we have 5(7) + /(/3) — 2 > 2(2 — 7) (5 — 7) (3 — f3). Thus the numerator of is 
positive. Since the denominator of is obviously negative, we get < 0. 

Hence since pm+i(2) ~ (—1)"'"''^ we have proved before, andpm+i(2)+Pm(2) ~ (— l)"^"*"^ 
by the induction assumption, we finally get ^^+2(2) + Pm+i{2) ~ (— 1)™. □ 
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Proof of Theorem A. Recall that in Lemma 4.4, we have proved that Pm+i {4>- (Am,r„)) ~ 
(_l)m+r,„-i g^^^ Pm+i{<P+i>^m,rJ) ~ (-1)"^+''-. Furthermore, A^+i^^^+i is the only root 

of Pm+l{x) between </)_(Am,r„i) &Ild </'+(Am,^r^). 

When m = 2, we have p3(0_(A2,r2)) > and ps (</)+( A2,r2)) < since r2 is odd. By 
Lemma A, we have ^3(2) > 0. Since X3,r2+i is the only root between </)_(A2,r2) ^^id 
0+(A2,r2), we get A3,^^2+i > 2. 

When m > 3, we have pm+i(0-(Am,r„)) ~ (-1)™"^ and prn+i(0+(Am,r™)) ~ (-l)"" 
since is always even. Still by Lemma A, we have Pm+i(2) ~ (—1)'""-'^. Since Xm+i^r„,+i 
is the only root between 0_(Am,r„) and 0+(Am,r„), we get A^+i^^^+i > 2. □ 

Remark. This theorem says that when m > 3, the first m-level initial eigenvalue is 
always greater than 2. 

Lemma B. Let m > 2. Then qm{x) > whenever < x < </)^^^(6). 

Proof. Define Omiz) = g^((/>L™^(z)) on < ^ < 6, Vm > 2. 

When m = 2, 62{z) = q2{(j)^^\z)) . Noticing that q2{x) = s{x) and q2{x) = —3x'^ + 

(2) 

22x — 35, an easy calculus shows that q2{x) is monotone decreasing when < x < 0i (6). 
Hence VO < x < 0L^^(6), we have ^2(0) = 26 > q2{x) > q2{<P^^\Q)) ~ 12.68. Thus 

26 > ^2(^) > 12.68,V0 < z < 6. (9.3) 

When m = 3, ^3(2) = ^3(0-^^)) = s{(f)^^\z))e2{z) - l{<j)'^^\z))r{(j)^^\z)) on < ^ < 6. 
Noticing that s{(f)^^\z)) = q2{4>-\z)) and q2{x) is monotone decreasing when < x < 
(P-\Q), we have 

s(0) = 26 > s{(l)^^\z)) > s((/)L^^(6)) ^ 22.96. 

(n\ (2^ 

The monotone property of — (2))r(0_ (2)) on < 2 < 6 implies that 

-84.21 > -l{(j)^^\z))r{(j)'-^\z)) > -120. 

Hence by using (19. 3p . we get 

26 ■ 26 - 84.21 = 591.80 > ^3(2) > 22.96 ■ 12.68 - 120 = 171.16, VO < z < 6. 

Hence ^3(2) > 6^2(^) > on < z < 6. 
We now use induction to prove: 

Om+i{z) > QOmi^) > on < 2 < 6, Vm > 2. (9.4) 

Of course, it holds for m = 2. To use the induction. Assuming 9m+i{z) > 69m{z) > 
on < 2 < 6, we will prove 9m+2{z) > 69m+i{z) > on < 2 < 6. 

Consider a polynomial 5'(x) = s(x) — |/(x)r(/(x)) = 6 + ^x — ^x'^ + ^x^ — ^^^ + l^^- 
It is easy to compute that 

g'{x) = - + 89x2 _ 64^3 ^ 5^4 > 115 _ ^^W^g) _ ^(0W(6))3 ^ 36.02 > 
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on < a; < C^(6). Hence g{x) is a monotone increasing function on the interval 
[0,0^^^ (6)]. So g{x) > g{0) = 6 on < x < 0L^^(6). 

By using an expansion along the last row of 6jn+2{z) = qm+2{4'^^~^'^\z)) , we have 

By the induction assumption and the fact that (f)^^'^'^\z) < (j)^I!^^^\z) < 2, we have 

= g{r^\z))e^+r{z). 

Since < 4>^^^'^\z) < 0^''(6) on < 2; < 6 when m > 2, we have g{(l)^I!^~^'^\z)) > 6. 
Hence 

em+2{z) > QOm+liz) > 0. 

Hence we have proved (19.41) holds for m + 1. From (19.41) . we get the desired result. □ 
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